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Abstract 

We call a state "vacuum bounded" if every measurement performed outside a specified 
interior region gives the same result as in the vacuum. Wc compute the maximum 
entropy of a vacuum-bounded state with a given energy for a one-dimensional model, 
with the aid of numerical calculations on a lattice. For large energies we show that 
a vacuum-bounded system with length Lin and a given energy has entropy no more 
than S'^^ + llnS'''^, where S'^^ is the entropy in a rigid box with the same size and 
energy. Assuming that the state resulting from the evaporation of a black hole is 
similar to a vacuum-bounded state, and that the similarity between vacuum-bounded 
and rigid box problems extends from 1 to 3 dimensions, we apply these results to 
the black hole information paradox. Under these assumptions we conclude that large 
amounts of information cannot be emitted in the final explosion of a black hole. 

We also consider vacuum-bounded states at very low energies and come to the 
surprising conclusion that the entropy of such a state can be much higher than that 
of a rigid box state with the same energy. For a fixed E we let L'in be the length of 
a rigid box which gives the same entropy as a vacuum-bounded state of length Lin. 
In the L' ^ limit we conjecture that the ratio L'in/ Lin grows without bound and 
support this conjecture with numerical computations. 
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Chapter 1 
Introduction 



1.1 Background 

Since the discovery of black hole radiation by Hawking [3], the fate of information 
faUing into a black hole has been a mystery. (See [4-6] for reviews.) If Hawking's 
semiclassical calculation is correct, then the outgoing radiation is purely thermal and 
the outgoing photons are uncorrelated to each other and to the matter which formed 
the black hole. If the evaporation is complete, and if the thermal nature of the 
radiation persists throughout the evaporation, then the original information is lost. 
Thus if the black hole is formed from a quantum-mechanically pure state, there will 
nevertheless be a mixed state after the evaporation is complete. This is the position 
held by Hawking (e.g. see [7]), but it violates CPT and may lead to difficulties with 
energy conservation and causality [4,8-10]. 

If information is not lost in black hole evaporation, there are several possibilities. 
One is that the black hole does not evaporate completely, but instead produces one or 
more Planck-scale remnants (e.g. see [6]). Another possibility is that the information 
disappears into a baby universe [11]. In this scenario the quantum-mechanical pure 
state is preserved, but parts of it are inaccessible to observation. It is also possible that 
the radiation is not really thermal, even at early times, because of a complementarity 
principle [12-15] or the inapplicability of the semiclassical approach [16-19], and thus 
that the information is encoded in subtle correlations in the radiation. In this case the 
black hole could act like a normal object with the entropy describing internal degrees 
of freedom. Some results from string theory [20-23] tend to confirm this view. 

Even if the radiation is thermal and uncorrelated during most of the evaporation, 
there is no reason to believe that it remains thermal near the endpoint of the evap- 
oration. The late-time radiation is presumably governed by an unknown theory of 
quantum gravity, and may well have correlations to the radiation emitted earlier.^ 
However, it is generally believed that late-time radiation cannot resolve the informa- 
tion paradox [24,5,6]. The argument goes as follows: While the black hole is large. 



^This is conceivable because at early times information in the outgoing radiation can be correlated 
with information in the ingoing negative-energy flux. 
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it is presumably radiating high-entropy thermal radiation. If the final explosion is 
to restore a pure state, it must radiate as much entropy^ as was radiated in earlier 
times. However, by the time the black hole reaches the point where unknown physics 
could come into play, there is little energy remaining. To radiate a lot of information 
with little energy requires a long period of time, and thus the "final explosion" looks 
more like a long-lived remnant. 

However, Wilczek and Holzhey [25] argue from a moving mirror model that a 
state with high entropy can nevertheless be purified with arbitrarily low energy cost. 
In certain ways, their model looks more like a remnant theory than a complete- 
evaporation theory, but it still appears to cast some doubt on the standard argument 
above. 

In any case, this argument requires bounding the entropy that can be contained 
in a particular region with a fixed energy. In the case of a region with refiecting walls, 
this is the question of finding the thermal state of quantum fields in a box. For a 
spherical box and a particular field theory the problem is easily solved. But with a 
region of complex shape, or where one wishes to make a statement intended to apply 
to all field theories, the situation is more comphcated. Bekenstein [26,27] argues that 
such a universal bound exists, but Unruh and Wald [28] disagree. 

Here we take a different approach. We consider only a single scalar field, but we 
use a weaker and, we hope, more physical condition on the results of the black hole 
evaporation. In the end, our results still support the claim that late-time radiation 
cannot restore the purity of the state of an evaporating black hole. 

We will also make a more general investigation of our new definition of a localized 
state. In the very low energy regime we will find the surprising result that this 
definition gives rise to much higher entropy than the same energy could support in 
a rigid box. In the low-energy limit we conjecture that this difference grows without 
bound. 

1.2 The vacuum-bounded state 
1.2.1 Black hole evaporation 

We start by considering a black hole formed from a pure quantum-mechanical state 
of incoming matter. To avoid any possible complications of quantum gravity theory, 
we will look at the state produced after the black hole has completely evaporated [5]. 
Gravity should play no significant role in this state, since the energy density should be 
small everywhere.'^ We can describe the final state as follows: at large distances from 
the position of the black hole (taken as the origin) there is outgoing Hawking radiation, 
which we are assuming to be thermal and without internal correlations. Within some 
distance R of the origin, there is some state of ordinary quantum fields that could 



^Here and throughout this thesis, "entropy" means fine-grained quantum- mechanical entropy. 
^If instead there are Planck-scale concentrations of energy, then we would have a remnant theory, 
a possibility we are explicitly not considering here. 
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have correlations with the radiation emitted earlier. The distance R is the distance 
that such information might have propagated since unknown physics came into play. 
Let us assume that Hawking's semi-classical calculation is good up to an energy scale 
Tunk. This temperature is reached when the black hole's mass is'' Mq = l/(87rTunk)- 
If, after this, the rate of evaporation continues to match the Hawking calculation,^ the 
black hole will evaporate in time t ~ IO'^Mq/q ~ l/{gT^^^), where g is the effective 
number of degrees of freedom in the particles that can be radiated. (See [29].) So 
there is a sphere of radius 



which contains total energy 



in which the information could be contained. Taking, for example, Tunk = lO^^GeV ~ 
lO^^mpi, and 51 ~ 100 we get^ 

R ~ 10^°/pi ~ lO^^^cm (1.3a) 
Eo ~ lO^mpi ~ lO^^g ~ lO^Vg. (1.3b) 

This yields a somewhat outrageous density of lO^^g/cm^, which is nevertheless small 
compared to the "GUT density" {lO^^GeVy/ (h^c^) - lO^^g/cm^ 

1.2.2 The vacuum-bounded condition 

Now we would like to answer the following question: How much entropy can be 
contained in a spherical region of radius R with a total energy EqI To answer this 
question we have to specify what we mean by "contained in a region." As mentioned 
earlier, if we ask how much entropy can be contained in a spherical box of radius 
R with perfectly reflecting walls, the question can be easily answered. However, 
the system with the box is not so closely akin to the system under discussion. For 
instance, inserting the reflecting walls into the system produces a divergent increase in 
the ground-state energy of the system. Furthermore, if we started with the vacuum in 
the whole system, and then introduced a spherical wall, we would produce a divergent 
geometric entropy [30-33]. A better description of our system is simply that it has 
thermal radiation outside radius R, and an unknown state of the quantum fields inside 
radius R, but no barrier or boundary at R. 

To study such systems, we will assume that the difference between the external 
Hawking radiation and an external vacuum is not important to considerations of 



^Wc arc working with units in which c=G = h=kB = i 

^As opposed, for example, to slowing to nothing and leaving a remnant. 

^Another possibility is that g diverges as T — > mp\. In this case the information can be radiated 
in a small number of particles of about the Planck mass, chosen from an infinite spectrum of such 
particles. This is effectively a remnant theory. 



R-z^ (1-1) 



unk 



Eo = ^ (1.2) 
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entropy. We will study systems that have an arbitrary state inside R but the vacuum 
outside R. To make this precise we will specify the problem as follows: 

Let a vacuum-bounded state be a generalized state (i.e. density matrix) 
for which every operator composed of fields at points outside a specified 
interior region has the same expectation value as in the vacuum. What is 
the maximum entropy of such a state whose interior region is a sphere of 
radius R and whose total average energy^ is given by (H) — EqI 

We expect to find that the answer to this question is similar to that of a box of 
radius R with reflecting walls, with some small correction. 

We will denote by X^^ quantities in the given interior region and by X^nt those 
outside this region. We will say that a generalized state is "localized to the inside" or 
"obeys the vacuum-bounded condition" if any measurement performed on the outside 
field operators in this state yields the same result as in the vacuum, i.e. if 



for every operator Oout that is constructed from field operators in the outside region. 

In the language of density matrices, we can write pout = TrinP, where p is the 
overall density matrix describing our system and Trin means to trace over all the 
"inside" variables. Then pout is the reduced density matrix describing only the outside 
variables, and Eq. (1.4) is equivalent to 



where |0) denotes the ground state. 



''If this approximation is bad we can increase R until the Hawking radiation outside R has very 
low temperature. 

^We cannot specify that every measurement of the energy must give Eq. Such a state is necessarily 
static and thus cannot represent outgoing radiation. 



TrpOout = Trp^Oout = (0|Oout|0) 



(1.4) 



pout = p:s^ = 'iVi„|o)(oi, 



(1.5) 



Chapter 2 

Mathematical Considerations 



In this chapter we will prove some theorems regarding the general theory of con- 
strained density operators. Our goals are to show that there is always a unique 
operator which maximizes the entropy subject to the constraints and to demonstrate 
the form of that operator. Unfortunately we will not accomplish either of these goals 
completely, although we will make some progress in these directions. In particular 
we will show that if such a operator exists it is unique, and that an operator of the 
form specified, if one exists, does in fact maximize the entropy. 

2.1 Topology of the space of density operators 

Definition Let H. he a separable Hilbert space, and let p be a linear operator on H. 
Then p is a density operator ( or density matrix) if and only if 

• p is Hermitian 

• p is trace- class and Trp = 1 

• {ip\p\ip) > for any state lip). 

The last property is usually written "p is positive semidefinite" (or sometimes "pos- 
itive definite") by physicists and "p is positive" or "p > 0" by mathematicians. In 
this section we will usually use the latter notation. 

We will require our operators p to satisfy a set of constraints of the form 

Tr pCa^Va (2.1) 

for some given operators Cq and numbers Va- (By TrA = a we mean that A is 
trace-class and its trace is a.) We would like to find the density operator p which 

maximizes the entropy S = — Trplnp subject to the constraints. There are two 
potential problems here. First, there might be states with arbitrarily large entropy. 
If there are no constraints, that will be the case whenever Ti, has infinite dimension. 
Second, there might be a supremum of the possible values of S which is not achieved 
by any p. 
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In an attempt to address these problems we would like to prove some continuity 
results for S and some compactness results for the space of allowable p. However, 
if any density operator is allowed, then those with infinite entropy or infinite energy 
are dense in the space of all density operators, using any reasonable norm. Thus S 
cannot be continuous on such a space. Instead we will work only with those p for 
which Tr pH is bounded, where H is some reasonable Hamiltonian. 

Since density operators are trace-class they are also Hilbert-Schmidt operators, 
and we can use the Hilbert-Schmidt inner product (p, p') — Tr pp' and the corre- 
sponding norm given by ||p||2 = ^/Tr~f^. The topology on the space of operators will 
always be that induced by ||-||2. We will write ||-|| for the norm of a state and for the 
bound on a bounded operator, and ||-||2 for the Hilbert-Schmidt norm of an operator. 

2.1.1 Compactness 

We will first show that with an energy bound the space of density operators is com- 
pact. 

Theorem 1 LefH he a separable Hilbert space, let H be a positive Hermitian operator 
on Ti with a purely discrete spectrum, and suppose further that no infinite- dimensional 
subspace exists on which H is bounded} For some E > let V{E) be the space of 
density operators p on Ti. such that TrpH < E. Then the space V is is compact in 
the topology induced by ||-||2- 

The proof is composed of a series of lemmas. 

Lemma 1 Let {An} be a sequence of operators on TC and suppose that Zim^-^oo^^"^ — 
A. Then MiUn^^A^^^ ~ Aa/s where A^p denotes a matrix element of A in some 
orthonormal basis {\cy)}. 

0. We can use the basis {|q;)} to compute 
^l/?)r = El^S-^«/5|'' (2-2) 

a/3 

Aap- I 

Lemma 2 Let {A^'^^} be a sequence of positive operators on Ti. with A^"") — > A. Then 
A>0. 

Proof: Let \x') be any element of ?i. Define a normalized vector \x) = |a;')/|||a;')||. 
We can always make a basis with |a;) as one of its elements. Prom lemma 1, A^ — > 
Axx- Thus 

{x\A\x) = lim(x|>l(")|x). (2.3) 



^If such an infinite-dimensional subspace does exist then there are states with infinite entropy for 
any E. 



Proof: If ^ A then -A^ 
the norm, 

|^H-A||2 = 5^|(a|^^") 

a/3 

SO each {Aap — A^^^l — >• and thus A^^^ - 



2.1. Topology of the space of density operators 
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Since A'-"^ > we have (xIA*^"-* |a;) > for any \x). Since the hmit of positive numbers 
cannot be negative, (xl^la;) > and thus (x'|A|a;') > for any That is to say, 
A>0.i 

Lemma 3 Let B be any positive operator on 7i with a purely discrete spectrum. Let 
{A*^"^} be a sequence of positive operators with A^"^ — > A. Suppose that TrA^'^^B < Bq 
for all n. Then TrAB < Bq. 

Proof: We work in a basis where B is diagonal. We write the trace 

N 

TtAB= hm y^A^A- (2.4) 

n=\ 

For each this is a finite sum. By lemma 1 the elements of AS^^ converge individually 
to the elements of A. Thus 

N N 

J2AnnBn=\imJ2A^lBn. (2.5) 

n=l n=l 

Since A^^^ > 0, the diagonal elements of A^'^^ cannot be negative, and since 5 > we 
have -B„ > for every n. Thus 

N 

Y,A^IB^<TtA'^^B<Bo. (2.6) 

n=l 

By Eq. (2.6) the terms in the limit in Eq. (2.5) are each no more than Bq, and thus 

N 

Y,KnBn<Bo. (2.7) 

n=l 

Thus each term in the limit in Eq. (2.4) is not more than Bq and so Ti AB < Bq. I 

Note that lemma 3 would not hold with < replaced by =, as shown by the following 
counterexample: Let B — diag(l, 2, 3, 4 . . .). For any n let A'^'"^ — except for the 
nth diagonal element which is 1/n. By making n large we can make ||A(")||2 — 1/n 
as small as desired. Thus A^'"^ A — 0. However, TrA^'^^B — 1 for all n, whereas 
Tr AB = 0. 

We will now work in a basis {|n)} where H is diagonal, Hmn = {m\H\n) = SmnEn 
with En < En+i. 

Lemma 4 Let A be any positive operator with TrAH < E. Given any e > there 
is some number N such that for every positive operator A with TrAH < E, 

oo 

^{n\A\n)<e. (2.8) 

n=N 
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Note that Tr AH < E implies that A is trace-class and consequently that for any 
given A, there is an N with J^'n "^nn < e- The point of lemma 4 is that such an N 
can be chosen uniformly for all A. 

Proof: We have 

oo 

J^AnnEn^TrAH <E. (2.9) 

n=l 

Since ^4 > we have Ann > for every n. Since the En are increasing, 

oo oo 

J2 AuuEn < Yl ^ ^ (2.10) 

n=N n=N 

SO 

°° _£! 
J2Ann<^. (2.11) 

Since H is unbounded on any infinite-dimensional subspace, we can find E^ arbitrarily 
large. Thus for any e > we can find an N such that 

oo 



J2^nn<e (2.12) 



n=N 



for every A with Tr AH <EA 

Lemma 5 Given any e > there is a 6 > such that for any positive operators A 
and A' with TrAH < E and TrAH' < E, if \\A - A'\\2 < 6 then \ Tr{A - A')\ < e. 

Proof: Let pIn denote the N x N matrix given by pab with a,h — 1...N. By 
lemma 4 we can find N such that 



Tr(^-^i^) = Y^Ann<{ (2.13) 



e 

n=N 

for any positive operator A with TrAH < E. Let 

S = . (2.14) 

Then for any A and A' with ||A — < 5 we have 

(AT \ 2 AT 2 

E(^- - < Nj^i^nn-A'nnf < N\\A-A'\\l < , 

n=l I a=\ 



(2.15) 



so 

\Tr{AU-NU)\<^-. (2.16) 
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If TrAH<E and TrA'H < E we can add Eqs. (2.13) and (2.16) to get 

\TrA-TrA'\ < |Tr (A - AU)| + |Tr (^ijv - A'iiv)| + |Tr (A'U - A')| < e.l (2.17) 

Lemma 6 Let A be a trace- class operator and C be a bounded operator. Then AC 
is trace-class and Tr\A\C < \ TrA\ x ||C||. 

Proof: Since A is trace-class it has a purely discrete spectrum. Thus we can 
write A = ^ A„|n)(n| with (n|n) = 1. By the triangle inequality, {n\C\n) < \\\n)\\ x 
||C|n)|| < ||C||. Then 

Tr\A\C^Yl \An\{n\C\n) < ^ (|^„| x ||C||) ^ \Tr A\ x \\C\\ . I (2.18) 

n n 

Lemma 7 Let {A^^^} be a sequence of positive trace-class operators with — > A 
and let C he any hounded operator. If TrA^'^^C — c for all i then TrAC — c. 

Compare with lemma 3. In the case of a bounded operator, taking the limit 
preserves the trace exactly, whereas with an unbounded operator we get only an 
upper bound. 

Proof: Since A^'^ A and using lemma 5, for any e > we can find \\A^''^ ~ ^lU 
sufficiently small that 



Tt{A^'^ -A)\ < 




(2.19) 



Then by lemma 6, Tr \A^^-A\C < e. Since Tr A^^C = c, it follows that |Tr AC-c\ < e 
for every e > and thus that Tr AC — c as desired. 

Lemma 8 The space V{E) is complete. 

Proof: Let {p*-*-*} be a Cauchy sequence of density operators in V{E). The space 
V{E) is part of the space of Hilbcrt-Schmidt operators on Ti, which is a complete 
space under the ||-||2 norm. Thus there is a Hilbert-Schmidt operator 

p = lim (2.20) 

i— >oo 

Clearly p is Hermitian. By lemma 2, p > 0. By lemma 3, Tr pH < E. Then by 
lemma 7 with C = /, Trp = 1. Thus p G V{E). So any Cauchy sequence in V{E) 
converges to a limit in V{E), which is to say that V{E) is complete. I 

Lemma 9 The space V{E) is totally hounded, i.e. any sequence of p^"'^ in V{E) has 
a Cauchy subsequence. 

Proof: Let {p(")}, n = 1.. . cxo be an infinite sequence of operators in V{E). As 
before we work in a basis where H is diagonal. The space of Hermitian, positive, unit- 
trace n X n matrices is compact. Thus any sequence of such matrices has a Cauchy 
subsequence. So define a set of sequences of integers {nk{i)}, k > as follows: Let 
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no{i) = i and for each k define {nk{i)} to be a subsequence of {nk-i{i)} such that 
jj. is a Cauchy sequence. Then let — We claim that p^'^ is a Cauchy 

sequence. 

Given any e > 0, by lemma 4 we can find an N such that 



oo 



(2.21) 



n=Ar 



for every p e V{E). Since p is positive, |pa6p 

oo oo 

EE 



IPabP < PnnP&6- ThuS 

oo oo oo oo 

Pab? < ^Pnn^Pbb - ^ ' ^ ^6 

a=l b=N a=l b=N 



(2.22) 



Now since {p^^^^'^Uat} is a Cauchy sequence we can choose a number M > N 
such that 



for all i,j > M. Then 



(2.23) 



E (njv(i)) _ (njv(i)) 
A'afe Pab 

a or b > N 

oo 00 

— |A^o6 Pab 

a=l b=N+l 



00 00 

E 



(0) 



,("iv(j)) 



Thus 



a=l b=N+l 

|p("Jv(«)) _ p{njv(j))||^ 



Pafc 



e 

< - 
2 



(2.24) 



IIP ~ ^ II2 — ^ * 
for all i,j > M. For any given k > N, nk{k) is n]y{i) for some i > k, since the 
are a subsequence of the n]s[{i). Thus for any k,l > N, and consequently fo 
k,l> M, we find that 



(2.25) 



(2.26) 



for some i, j > M, and thus that 

llpC^) -p«||2 <e (2.27) 

which is to say that p^*^ is a Cauchy sequence. I 

Proof of theorem 1: By lemma 9, V{E) is totally bounded. By lemma 8, V{E) 
is complete, which is to say that V{E) is compact. I 



2.1. Topology of the space of density operators 



23 



2.1.2 Continuity of the entropy 

We now prove that S is continuous if we restrict ourselves to p with Tr pH < E. 
Theorem 2 S{p) — — TVplnp is continuous on V{E) in the ||-||2 norm. 

Proof: Let p be any density operator in V{E). Given any e > we need to find 
a 5 > such that \S{p) - S{p')\ < e for every p' G V{E) with \\p - p'y < 6. We will 
do this by writing p as the limit of N x N matrices pjjv and p' as the limit of p'In- 

We work again in a basis where H is diagonal. Since pl^ is a matrix of finite 
dimension, S{pIn) is continuous. Thus given any N we can find a S such that 



whenever WpIn — p'InIU < ^- Thus it is sufficient to find an N such that \S{pIn) — 
S{p)\ < e/3 for all p G V{E). 

Let Pk{p) denote the kth largest eigenvalue of p. Lieb, Ruskai and Simon [34] 
show that Pk{piN) < Pk{p), and go on to show from this that 



We want to extend this result to show that S{pl]y) — > S{p) uniformly for all p such 
that Tr pH < E. 

Define the function 



Note that s is continuous. It is increasing for p < 1/e and decreasing for p > 1/e. As 
we go from S{plj\f) to S{p) the pk will increase. It is possible to have a few large pk 
such that increasing p^ will decrease s{pk), but for the majority of an increase in 
Pk means an increase in s{pk). 

First we study how much S can decrease as N —>■ oo. Let e be any small positive 
number. By lemma 4 we can find N_ large enough so that 



\s{piN)-s{p'u)\<'- 



(2.28) 



hm S{pU) = S{p) . 



(2.29) 




(2.30) 



oo 




(2.31) 



a=Ar_ 



for every p eV{E). Thus 



oo 



Tr pIn = 1 - ^ Pnn > 1 



e 



(2.32) 



3 



n=N+l 



for any N > N_. Now ^Pk{p) = Trp = 1, so 



^{Pk{p) -Pfe(pijv)) < o ■ 



(2.33) 
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Since ds/dp> —1 for all p < 1, we can conclude that 

S{piN) - S{p) < 

for all N > N^. 

Now we study how much S can increase as N ^ oo. We write 



(2.34) 



P 



pIn 


A 







(2.35) 



For fixed p[n and pt , S{p) will be largest when A = Q. To see this, consider 



pIn 


-A 


-A^ 





(2.36) 



The eigenvalues of p_ are the same as those of p, so S'(p) — S{p_). Now let 



Po 



pIn 


^ 







pr J 





(2.37) 



Since the entropy is convex (see [35] and section 2.2.3), S{po) > {S{p) + S{p^))/2 
S{p). Thus S{p) < S{po) = SipU) + -S(pT^). 
Now pt^ obeys the constraints 



TVpr = l-TVpU = l-p< 



(2.38a) 
(2.38b) 



Let Sth{E; ei, 62, 63, . . .) denote the thermal entropy of a system whose energy levels 

are ei, 62, 63, Suppose we had a system with energy levels 0, E'at+i, -E'Ar+2, ■ ■ ■ ■ 

Then a possible normalized density matrix would be 



(2.39) 



p< 








pr 



It would have energy E^ < E and entropy — p<lnp< + S{p^^). This must be less 
than the maximum entropy for this system with energy < E, so we conclude that 



S{p) - S{pU) < 5(pT^) < ^th(^; 0, ^iv+i, ...). 



(2.40) 



Now as A?" — > 00, only higher and higher energy states are permitted in the expression 
above, so the entropy will decrease to zero. That is. 



lim Sth{E; 0, En+i, ^7v+2, . . .) = . 

JV— »oo 



(2.41) 
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So, given e > we can find A'^-i. such that if > then 



S{p) - SipU) < 



(2.42) 



for any p eV{E). 

Now we set N = max(iV+, iV_). From Eqs. (2.34) and (2.42) we find that 



for every p e V{E). Then we choose (5 as in Eq. (2.28). If \\p — p'\\2 < S then 



\S{p)-S{p')\ < \S{p)-S{piM)\ + \S{piN)-S{p'iM)\ + \S{p'iN)-S{p')\<e.i (2.44) 



Theorem 3 Let Ti. and H he defined as in theorem 1, let {Cq} he a set of positive 
operators with discrete spectra, and let V<{E; Vi, V2, . . .) he the set of density operators 
p such that TrpH < E and TrpC^ < for all a. Then S{p) achieves a maximum 



Proof: From theorem 1 we know that V{E) is compact. Thus any sequence of 
pW £ V^(E;Vi,V2, . . .) has a subsequence that converges to some p e V{E). By 
repeated apphcation of lemma 3 we see that TrpC^ < so p e V<{E] Vi, 1^2, ■ ■ ■)■ 
Thus V<{E; Vi, V2, . . .) is sequentially compact, and so it is compact. By theorem 2, 
S is continuous on V{E)., so S achieves a maximum on 'P<{E\ Vi, V2, . . .). I 

It is possible to rearrange the constraints that we will need into positive operators. 
However, our constraints are of the type Tr pCa — Va rather than Tr pCa < Va- What 
we really need is 

Conjecture 1 Define Ti, H and {Cq} as in theorem 3 and let V={E; Vi, V2, . . .) he 
the set of density operators p such that TrpH — E and TrpCa — for all a. Then 
S{p) achieves a maximum on V={E; Vi, V2, . . .). 

Unfortunately, we do not know how to complete the proof of this conjecture. However, 
in the next section we show that if there is such a p it is unique. Furthermore since 
V={E; Vi, V2, . . .) C V<(E\ Vi, V2, . . .) the maximum value of S in theorem 3 is at 
worst an upper bound on S{p) with p e V={E; Vi, V2, . . .). 



To maximize the entropy we would like to look at the variation of 5 as p is varied and 
require that ^yS" = so that S can be a maximum. However, when p is varied we may 
find that SS is not well-defined. There are two possible causes for this. One is that 
p has some zero eigenvalues. The second is that p has arbitrarily small eigenvalues 
which become negative under any variation. We treat these problems below. 

First, since we are concerned with zero eigenvalues we should specify what we 
mean by S. When we write S{p) = — Tr plnp we mean that S{p) = s{p), where s is 



\S(p)-S(pU)\<^, 



(2.43) 



onV<{E;V,,V2,...). 



2.2 Maximization of the entropy 
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defined in Eq. (2.30). That is, if p = '^Pa\a){a\ tlien S = ^s(Pq) or alternatively 
if we write s{p) as a power series in p then S{p) is given by the same power series in 
p. Under this definition S{p) is well-defined when p has zero eigenvalues, even though 
In p is not defined in such cases. 



We want to vary p via p' = p + t5p and look at the resulting change in S{p'). We will 
call a variation Sp admissible if: 

• 5p is Hermitian. 

• Tr(5p = 0, so that Tr p' = 1. 

• For each a, Tr 5pCa = 0, so that Tr p'Ca is unchanged. 

• For sufficiently small t > 0, p' > 0. 

Under these conditions p' is a well-defined density operator for at least some small 
range of t > 0. 

We will often want to let Sp interpolate between p and some other density operator 
p. To show that this is always possible, we prove 

Theorem 4 The space of density operators p that obey the constraints is convex, 
i.e. if p and p are density operators that obey the constraints then Sp = p — p is an 
admissible variation and po — p + tSp is a density operator for all t e [0, 1] . 

Proof: 

• Since p and p are Hermitian, so is Sp. 

• Tr5p = Trp-Trp = 0. 

• Tr SpCa = Tr pCa - TVpC« = 0. 

• Let X be any state vector. Then for all t E [0, 1], x • p'x — x - {tp + {1 — t)p)x = 
t{x ■ px) + (1 - t){x ■ px) > 0, i.e. p' > 0. 1 



where pa are the eigenvalues of p. When we vary p, the change in the eigenvalues is 
given by first-order perturbation theory. 



2.2.1 Variation of S 



Now 




(2.45) 



a 



dp, 



'a 



{a\Sp\a) = Sp, 



(2.46) 



dt 



aa 



If Pa > then we have 



dsjpa) 
dpa 



(l + lnpj. 



(2.47) 
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Thus if there are no zero eigenvalues we can write 
dS 

— = - + lnp«)(^Pa = -Tr 5p - Tr 5plnp' = -Tr Spin p' , (2.48) 

a 

since we have required TrSp — 0. 

2.2.2 Prohibited sectors 

If there are zero eigenvalues, then Inp and thus Eq. (2.48) are not well-defined. We 
handle this case as follows: 

Theorem 5 If p maximizes S subject to the constraints TrpCa = and if p\oc) = 0, 
then every p that satisfies the constraints must annihilate \a) . 

Proof: Suppose to the contrary that p\a) = but there is some p such that 
Tr pC a = Va but p\a) ^ 0. We let 6p = p — p. By theorem 4, 6p is admissible. Now 
we look at just those Pa that are in fact changed by 6p. Let t be some small but 
positive number. Then all the Pa that are changed by 6p are positive and we can 
write 

^ = - E '^^^l"^- (2-49) 
a with Spa 7^ 

As t — > the — Inpo: terms grow without bound for those a where Pa{t) 0. Thus 
eventually these terms dominate everything else in Eq. (2.49) so that 

lim — = +00 . (2.50) 

t^o dt 

Thus for sufficiently small t > there is a p with S{p) > S{p) in contradiction to 
assumption. Thus if p maximizes S then we must have p\a) =0 for every \a) where 
p\a) = 0. 1 

Thus p can have a nontrivial nuUspace A/" C 7i only if every every p that satisfies 
the constraints annihilates the space A/". In this case, we will write 7i = A/" © and 
work in the space 7i'. We will replace all our operators with ones that act only on 
7i'. If, restricted to H', the constraints are linearly dependent, we can now discard 
some of them to have again a linearly independent set. This process eliminates the 
troublesome sector from the problem. 

2.2.3 Derivatives of S 

If Sp is an admissible variation, and as long as there are no zero eigenvalues, we can 
differentiate S with respect to t as in Eq. (2.48). We demand that 

^ = -Tr(5plnp' = (2.51) 
dt 
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for 5* to be a maximum. 

For the second derivative we care only about the sign. Lieb [35] showed that S is 
always (downward) concave. Here we obtain the same result by a different technique 
and show also that the concavity is strict, i.e. that 

§<0 (2.52) 

for any p and any nonzero Sp. 
Prom Eq. (2.51) 

To expand this we use the formula 



I In A = j\s {I- s{A s{A - I))-' 

which can easily be derived as follows: Let B — A — I. We expand 



and differentiate to get 



We observe that this is 



as desired. 

Using Eq. (2.54) we find 



dt^ 



(2.54) 



B^ 

lnA = ln(7 + S) = S + — + — + ••• (2.55) 



d , , dB 1 fdB ^ dB 
— ln^= -|- + - —B + B— 
dt dt 2\ dt dt 

+ -(—B^ + B—B + B^—] + ---. 2.56 
3\ dt dt dt ' ^ ' 



\ dB dB ^ dB\ 2 dB dB^ r.2dB\ 

ds— + s —B + B— + —B^ + B—B + B^ — 
dt \dt dt J V dt dt dt J 

= /' ds {1 + SB + s^B^ + ...)^(i + sB + s^B^ + ■■■) 
Jo dt 

= r ds {I - sB)-' ^ {I - sB)-' (2.57) 
Jo dt 



- [ dsTr5p{I - s{p - I))-^5p{I - s{p - I)y^ 
Jo 

- I TrX{sf (2.58) 

^0 



2.3. The form of p 
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where 

Xis) ^ {5pf'\l - sip - I))-\5pY'' . (2.59) 

Since X[s) is Hermitian, TrX(s)^ > with equality obtained only for X[s) = 0. 
Thus 

§ < (2.60) 

with equality only if X{0) = for all s. But -^(0) = Sp, and thus 

^"^^ < (2.61) 



for any nonzero 5p. 

Now suppose that p maximizes S. Then for any admissible 5p, we must have 

or else for sufficiently small t we would have S{p') > S{p). Furthermore we can see 
that condition (2.62) is sufficient to show that p maximizes S, as follows: Let p be 
any other density operator that meets the constraints. Let Sp = p — p. From theorem 
4, 6p is admissible. From Eqs. (2.62) and (2.52), it follows that S{p) < S{p). Since 
this holds for any p, p is the global maximum of S. It follows from this that if there 
is a p in conjecture 1 it is unique. 

2.3 The form of p 

We must find the unique state p (if any) which satisfies our constraints and which 
gives TrSplnp = for any Sp which maintains the constraints. This means that we 
are concerned with Sp such that 

Tr 5p = (2.63a) 
Tr SpCa^O (2.63b) 

for all a. We have included the Hamiltonian among the constraints Ca- For every Sp 
satisfying Eqs. (2.63) we must have 

TiSplnp^O. (2.64) 

If we choose 

In p = const + ^ faCa (2.65) 

a 

which is to say 

poce^f"^", (2.66) 
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with any coefficients fa, then Eqs. (2.63) ensure that Eq. (2.64) is satisfied. If we can 
find a p with the form of Eq. (2.66) that satisfies the constraints, then we have found 
the unique solution. We have not shown that such a solution always exists, but in 
the numerical work described in section 4.3 we have always succeeded in finding one. 
Since Tr p = 1, we can write 

P = lVeE/«c„ • (2-67) 

Our goal is now to determine the coefficients so that the constraints are satisfied. 
We can define a grand partition function, 

g^T^gE^Ca_ (^2.68) 

Its derivatives are 

^ = TrC„eS^"^«, (2.69) 

so 

(C«) = -^lng. (2.70) 
We have the usual thermodynamic formula for the entropy, 

S = -(Inp) = InQ - ^ UCa) . (2.71) 
Differentiating this we find 



Now we specialize to the case where one of the constraints is just the Hamiltonian. 
The corresponding coefficient is written — /3, and we have 

p^l^^-PH+Y.UC^_ (2.73) 

If we vary the coefficients in such a way that {H) = E changes but the other {Ca) = Vo. 
remain fixed we see that 

dS^-Y,jad{Ga)^m- (2.74) 



Thus /5 = dS/dE and so the coefficient /3 has the usual interpretation as the inverse 
temperature. 



Chapter 3 
The Problem 



Now we return to the problem at hand. We would like to find the maximum-entropy 
vacuum-bounded state with a given average energy Eq. We first reduce the problem to 
a simpler case which we hope will capture the important behavior. Then we analyze 
the sectors on which p must be zero and introduce new coordinates in which these 
sectors do not appear. We work with these new coordinates to derive formulas for 
the expectation values of the operators that must agree with the vacuum. Finally we 
look at the form of the results to be expected when we solve the problem numerically 
in the next chapter. 

3.1 Simplifications 

3.1.1 One scalar field, one dimension 

First we restrict ourselves to a theory consisting only of gravity and one massless 
scalar field. In such a system we can imagine preparing an incoming shell of scalar 
particles to form a black hole, which would then evaporate by emitting scalar quanta. 
Thus we can ask the same questions about black hole evaporation in this system as 
we could, say, in the standard model plus gravity. 

We will also begin here by working in one dimension. We will put our entire 
system in a box of length L and require that all deviations from the vacuum are in 
the region from to Lin. Later wc will take the overall box size L to infinity while L[n 
remains fixed. The inside region will be [0, Lin] and the outside region will be [Lin, L]. 
We will use the usual scalar field Hamiltonian, which in classical form is 
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3.1.2 Gaussian form 

In our problem, the constraints are the energy bound, 

TrpH^Eo, (3.2) 

and the vacuum-bounded condition, 

IVpO:,, = (0|O:,,|0), (3.3) 

where O^^^ is any operator which is constructed out of the fields 4>{x) and 7r{x) in the 
outside region. We will assume the solution has the form 

poce-^^+^^"0°"t. (3.4) 

We now show that fa is nonzero only for those operators O^^^. which are quadratic in 
the fields. 

Suppose that we wanted to solve a different problem in which we cared only about 
the constraints involving the quadratic operators. We would have the energy bound 
and the constraints 

Tt pQ:,, = {0\Q: JO) , (3.5) 

where Qout runs only over quadratic operators (f){x)(f){y) and 7!'{x)ir{y). (The operators 
(f){x)7i{y) vanish automatically by symmetry under —>■ —cj).) We expect the solution 
to this problem to have the form 

p'oce-^^+S-^-^^-t. (3.6) 

Now p' is a Gaussian operator; i.e., (<;/'(•) |p'l0'(")) ^ Gaussian functional of the 
values of (j) and (j)' . Let p'^^^ = Tti^ p'. The trace is just a set of Gaussian integrals, 
which means that the resulting p',^^ is also a Gaussian. Because H is quadratic, 
the vacuum p™'^ = |0)(0| is Gaussian, and so is its trace p™^ = TrinP™'^. Now 
by construction we have Tr p'^^^(f){x)(l){y) = Trp^0(x)0(y) and Tr p'^^^7r{x)7r{y) = 
Tr pl^^TT{x)7:{y). These conditions are sufficient to fix the coefficients in the Gaussian 
Po^^, and thus to show that in fact p^^^ and p^^ are the same Gaussian; i.e. that 
Trinp' = TrinP™^ 

Thus p' satisfies all the constraints of the original problem. Since only one p can 
have these properties it follows that p — p' and thus that Gaussian solution p' is the 
correct solution to the original problem. 

3.1.3 The discrete case 

We now approximate the continuum by a one-dimensional lattice of coupled oscilla- 
tors, with a classical Hamiltonian 



// = i(P,-P, + x-i^x) 



(3.7) 



3.2. Prohibited sectors 
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The simple kinetic term in Eq. (3.7) corresponds to choosing oscillators of unit mass, 
regardless of how densely they are packed. In terms of scalar-field variables this means 
that X/^ — ^/Ll(f)^ and = tt^/ \/Ll where Li — L/ {N + 1) is the lattice spacing, 0^ 
is the average of over an interval of length Li, and tt^ is the total momentum 
7r(x) in the interval. 

The matrix K gives the couplings between the oscillators and represents the d4>/dx 
term in the scalar field Hamiltonian. To approximate the continuum with the zero- 
field boundary condition we will imagine that we have N oscillators located at the 
points l/(iV + 1) . . . A^/(A^ -I- 1) and that the end oscillators are coupled to fixed-zero 
oscillators at and 1. Then 



K ^ 








-9 













v 



(3.8) 



where g = l/Lf. 

We will take A^in ~ i^in/^i of the oscillators to represent the inside region, and 
-^out — N — to represent the outside region. 

We want to maximize S subject to the constraints 



Tr pH = Eq 

Ti pXiXj = {0\xiXj\0) 
TvpPiPj = (0|PiP,|0), 



(3.9a) 
(3.9b) 
(3.9c) 



where i and j run over the oscillators which represent the outside region. We will 
define matrices X and P whose elements are the quadratic operators via 



X/Ui/ — 



P = 



so that Eqs. (3.9b) and (3.9c) become 

Tr pXout,out 

,out 



p p 



(0|Xo,t,out|0) 
(0|Pout,out|0) . 



(3.10a) 
(3.10b) 



(3.11a) 
(3.11b) 



3.2 Prohibited sectors 

As discussed in section 2.2.2, the first thing we need to do is to look for sectors of our 
Hilbert space that any p must annihilate in order to meet the constraints, and restrict 
our attention to the subspace orthogonal to these. We will start by examining the 
ground state of our system and looking for sectors which are forced to remain in the 
ground state. 
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3.2.1 A different description of the vacuum 

To work with the Hamiltonian of Eq. (3.7) we will make a change of coordinate to 
put it in diagonal form. Let Z he a. matrix whose columns are the eigenvectors of K, 
K — ZQqZ~^, with the normalization 

ZQoZ'^^I and Z'^^ftoZ'^ ^ K . (3.12) 

Define new coordinates z via x = Zz and P^^ = P^. In these coordinates, 

h=Ij:^':hpi+zI). (3.13) 

a 

The vacuum is the ground state of this Hamiltonian. We can define raising and 
lowering operators 

aa = {za + iPzJ (3.14a) 

ai = ^{^a-tP.J (3.14b) 

H-j:^^:'Ua + l) . (3.14c) 

The vacuum is the state |0) annihilated by all the a^. It is straightforward to write 
the expectation values in the vacuum state, 

= {0\zaZf3\0) = ^S^p (3.15a) 

{F,)ap^{0\P.^P.,\0) = ^Sap (3.15b) 

so 

(0|X|0) = Z{0\Z\0)Z^ = -ZZ^ (3.16a) 

(0|P|0) = Z-'^{0\F,\0)Z-^ = ^Z-^^Z-' = ^((0|X|0))-^ . (3.16b) 

There can be many different Hamiltonians that have the same ground state. If 
we consider 

//'=^(P,-T'P, + x-ir'x) (3.17) 

with T' and K' some coupling matrices, we can follow the above derivation to get a 
normal mode matrix Y and some frequencies Q with 



YQY'^ = T' (3.18a) 
y-^^Qy-^ = K' (3.18b) 
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X = Fy (3.18c) 
= Y-'^Py (3.18d) 

H'-lY.^p{P^+yp)- (3-18e) 



We can define raising and lowering operators for these modes, 

1 
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'p^^{yp + tPy,) (3.19a) 



bl = ^{y(^-^Py,) (3.19b) 
H'^J2^p(b^b + ^^ . (3.19c) 

Now y = Y~^K = Y~^Zz = W~^7. where 

W = Z-^Y . (3.20) 
Similarly, P^, = Y^P^ = Y^Z-^^P, = W^P,. Consequently 

= I {W^ai^a + ai) + Wlia^ - ai)) 

= I {{W-' + W^)paaa + - W^)f3aai) . (3.21) 

For H and H' to have the same vacuum we require that the bp depend only on the 

and not on the a|^, which is to say that — W'^, i.e. that is a unitary matrix. 
With W unitary YY'^ = ZZ^ so the vacuum expectation values of Eq. (3.16) have 
the same values expressed in terms of Y as they had in terms of Z. 

3.2.2 Modes that remain in the ground state 

Now consider a unitary matrix W and let Y = ZW as in the last section. Suppose 
we can find W such that Y has the following property: 

The A'" modes can be divided into A^'gs > "ground state" modes and 
A^frcc = N—Ngs "free" modes such that for all a < Ni^ and for all > A^free, 
Yaf3 = and Y^-' = 0. 
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That is to say Y and Y ^ will have the form 



Y 







^n,free 





^ut,free 


^ut,gs 




(3.22) 



and 



Y- 







Nout 






free, in 


free, out 






V " 


gs,out 





(3.23) 



If this is the case, then for (3 > A^free; bp — {Yp^x^ + lY^pP^^^ /y/2 depends only 
on outside operators Xi and P^^.. In any vacuum-bounded state, regardless of entropy 
considerations, (xiXj) and {PiPj) have their vacuum values. Consequently, for all 
(3 > A^free, if P describes a vacuum-bounded state then Trp6^6^ = (0|6^6/3|0) = 0. 
The vacuum-bounded constraint forces modes /3 > A^f^e to be in their ground states. 
These modes will not contribute to the calculation of the maximum-entropy vacuum- 
bounded state. 

How many such modes can exist? Let Wp denote a column of W with > Nf^ee 
and let Z"- denote a row of Z with a < iVin. Similarly let {Z~^)a denote a column 
of Z~^. Since Y — ZW, Yap — whenever Wp is orthogonal to Z". Similarly 
Y~^ = W^Z~^ so Yg~^ = whenever Wp is orthogonal to {Z~^)a- Since W is unitary, 
the Wp must also be orthogonal to each other. Thus there are iVgs columns of W 
which have to be orthogonal to A^m rows of Z, to A^in columns of Z~^, and to each 
other. Since there are A^ components in a column of W it can be orthogonal in general 
to at most N — 1 other vectors. Thus A^gs is limited by A?" — 1 = 2A^in + A^gs — 1, or 
= N-2Ni^ = Nont-Nin. Thus whenever A^out > there will be A^gs = A^out-A^in 
modes that are forced to remain in the ground state. 

These conditions determine the columns Wp with f3 > 2Ni^-, up to a unitary trans- 
formation on these columns alone, and likewise the remaining columns are determined 
up to a unitary matrix which combines them. 



3.2.3 Density matrix and entropy 

We can write Ti = Ti ® Q where is the Hilbert space of states of the "free" modes 
and Q is the Hilbert space of states of the "ground state" modes. Let |0)gs G Q denote 
the ground state of this system. For any operator A we can define an operator A that 

acts on ii via {c(\A\(3) = {a® Ogs|A|/? ® Ogs) for all a,(3 eH. 

Let p describe a vacuum-bounded state and let pgs = Tr frceP- Then Tr p^h^^hp — 0. 
This defines the vacuum state, so pgs = |0)gs(0|gs and thus 



P = P(g)(|0)gs(0|gs) 



(3.24) 
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If we write p in diagonal form, 

p = ^P,|<5)((5|, (3.25) 

a 

then 

a 

The entropy is 

,5 = -Trplnp= -^s(P„) = -Trplnp. (3.27) 

a 

Since any vacuum-bounded p has this form, a variation of p that preserves the con- 
straints must also have this form, 

= (|0)g,(0|gs) . (3.28) 
If ,5 = -Tr (p + tSp) ln(p + tSp) = -Tr {p + tSp) ln(p + tSp), then 

— = -Tropin p. (3.29) 

If p maximizes S subject to the constraints, then we must have Tr^plnp for any 
variation 5p that preserves the constraints, i.e. for which Tr5p = and 

Tr SpH = Tr SpXiXj = Tr SpPiPj = . (3.30) 

Now for any A, Ti 6pA = Tr (5p ® (|0gs)(0gs|))A. When we take the trace we only 
need to sum over states of the form \a (8) Ogg). Thus 

(A) = Tr(5pA= ^(a|(5p|;5)(/3(g)0gs|A|a(8)0gs) = Tr(5pA= (i). (3.31) 

dp 

Thus wc a£e looking fo£ p_that maximizes 5* = — Tr pin p subject to the constraints 
Tr 6pH = Tr 6px^j = Tr SpPiPj = 0, where i and j range over the outside oscillators. 
There is no longer a problem of zero eigenvalues of p. 

From Eq. (2.66) we expect p to have the form 

p (X e-f^i^+fij^i+aijPiPj) _ (3.32) 
We can write this in a more familiar way as 

p (X e-^"' , (3.33) 
where H' is a fictitious Hamiltonian for these oscillators, 

H'^H + fij^j + QijP^j . (3.34) 
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3.2.4 New coordinates 

We would now like to introduce new coordinates w as follows: The first A^m w coor- 
dinates will be the inside oscillator coordinates, Win = ^in- The last A^gg coordinates 
are the ground state normal modes, Wgg = yg^. The remaining N^^ coordinates can 
be any coordinates that are independent of those specified so far; we will make a 
particular choice later. 

To do this, we proceed as follows: From Eqs. (3.22) and (3.23) we have 



yy-l 





^in,free^free,in 


^in.frec^frccout 




i ^ut,free^ree,in 


^ut,free^ree,out "1" ^out,gs^s,out 



and 



-'free.in-^m.frec "r frcc.out out, free 



KS.out -' out, free 



V-1 V 
-'free.out-'out.gs 



gs,out-' out,gs 





iVout 









/ 


-^free 















}jVin 

(3.35) 



(3.36) 

In particular, l^g7,out^out,gs = I- We would like to extend 1^7,out ^ut,gs into square 
matrices R and with 



and 



{ D \ yout,gs )}jVout 

}Afout 



R~' 



iVout 

D' 

gSjOUt 



This means that we must find D and D' such that 

-'gs,out-^ ~ 
-D'l'out.gs = 

D'D = 7 . 



(3.37) 



(3.38) 



(3.39a) 
(3.39b) 
(3.39c) 



There are many possible choices of D and D' that satisfy Eqs. (3.39). Here we proceed 
as follows: Let 



Z = 



Zout 



(3.40a) 
(3.40b) 
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so that ^in^in^ = I, ^out^out = ^ ^nd Z-^^Zin + ^out^out = I- Now let 

D = ^Zout^S - (0|X|0)out,in (3.41a) 

D' = ^Zr^^Z„-i = (0|P|0)in,out , (3.41b) 

let A and B be A^in x Ni^ matrices with 

AB = (D'D)-^ (3.42) 

and let 

D = DA (3.43a) 

D' = . (3.43b) 

We also divide 

W = ( M^free | W^gs ) . (3.44) 

Since Y ^ ZW and = W^Z-\ we have 



ZinWgs = (3.45a) 

gs m 



W^,Zr^ = . (3.45b) 



Thus 



D%.t,,s « 5Zri^Z„-iZ,,tl^g3 = BZrfw,, - BZrf Zr^'Z,,W,, = (3.46) 

by Eqs. (3.45) and their transposes. Similarly 

« = PFjZ.^A - PFjZr^Zi^CA = . (3.47) 

From Eq. (3.42) wc find D'D = I. Thus the matrices D and D' satisfy Eqs. (3.39). 
We still have the freedom of choosing the matrix A arbitrarily. 
Now let 



Q 





iVout ^ 














R 


1 }Nout 



(3.48) 



and define w by x = Qw so P^; = P^. Then = Xjn and Wgg = ygg as desired. 
3.2.5 Reduced operators 

We would like to recast our problem in terms of Wfrcc, the the first A^frcc w coordinates. 
First we look at the operators x^x^, = X^^, and P^Pu = ^^u- If we write 

Q = (<5free|Qgs) (3.49a) 
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(3.49b) 



we have X = Qw = QfreeWfree+QgsWgs and P = Q ^ = Q^^l^ {Pw)fTee + Qgs(Pw)gs, 

so 

X = (3freeWfree,freeQfree + QgsWgs,gsQgs (3.50a) 
P = Qt£ (P..)free,freeQfTee + Q^f i^w),s,,.Q-s ■ (3.50b) 

Now we form the reduced operators and PfiPu. Since Wgs = ygg and the ygg 
modes are in the ground state by definition, Wgg^gs = (IPM!)gsgs = (1/2)/, and so 

X = gfreeWfree,free<5L + ^QgsQgs (3-51a) 
P - gSe^ (F.)free,freegf-e'e + ^Qg"f Qgs' " (3-51b) 

In each case the second term is just a constant. Now 

H^^ {P^P^ + K^^x^x,) = ^IV (P + KX) (3.52) 
where the trace is over the oscillator indices. Thus 

H = ^Tr (^(3f^eJ(P^)free,freeQfrle + KQhee^free,ireeQlee) + COUSt 

= ^Tr (T'(P^)free,free + ^Wfree.free) + COUSt (3.53) 

where 



K = Qlree^Qfree = 

















( 


K- ■ 



out, in 



-^in,out-D 



The constant term in Eq. (3.53) is 

-Tr ^KQgsQ^^ + ^ = -Tr ^i^ut,gs-^out,out^ut,gs + ^s,out^gs,out) 



(3.54a) 



(3.54b) 



(3.55) 



It depends on the ground state modes only and is just part of the zero-point energy. 
It will be the same in the vacuum and in a vacuum-bounded state. Thus if instead 
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of Eq. (3.53) we use 



-Tr 
2 



T(P^)free,free + ™: 



free, free 



(3.56) 



we are just shifting if by a constant term and thus changing the zero-point energy. 

Now we would hke to make this reduced system look as much as possible like the 
system we started with. Let Bi be some matrix such that BfBi = {D'D'^)^^, let 
i?2 be a unitary matrix to be determined, and let B = B2B1. Then D'D'^ = I, so 



T = 7. This gives A 
and let 



Now 



741^42 where Ai 



{D'D)-^B-^ and A2 



in,out 



-f^out,out-Dl 







Bj. Let Di = DAi 



(3.57) 



-9\ 



(3.58) 



so i^in,out-Di is nonzero only in the last row. Thus the last iVin + 1 rows and columns 
of Ki look like 



1 










1}' 






}iVi„ 



A matrix of size (iVin + 1) x (iVji 
transform of the form 



(3.59) 



1) can be put in tridiagonal form by a unitary 





^ 


1}' 









(3.60) 



which can be constructed using the Householder process. (See, e.g., [36] section 11.2.) 
We will use this to choose A^ — U, so that K will be tridiagonal. For each off-diagonal 
element in the resulting tridiagonal matrix there is a choice of sign, and we will choose 
them all to be negative. Thus in the w coordinates, each oscillator has unit mass and 
is coupled only to its neighbors. This completely specifies the matrices A and B and 
thus D and D' . Note that D and D' do not depend on the undetermined parts of W. 

To make H' in Eq. (3.34) we can add to ^ a kinetic and potential term involving 
outside oscillators only. The potential term (disregarding a constant) is 



fijX^Xj Tr /Xout,out Tr /QquI; free ^Wf].ee,freeQout, free 



Since 



Qout,free = ( | -D ) }jVout , 



(3.61) 



(3.62) 



this term is equivalent to adding an arbitrary term to just the lower right part of K. 
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Similarly the kinetic term is 



QijPiPj — Trfi'Pout.out — Tr5'Qfree,out(^ty)free,free<5free 



out 



Since 



QfrL 



iVout 

\ jiVout 



free,out 



(3.63) 



(3.64) 



this term corresponds to adding an arbitrary term to just the lower right part of f. 
That is to say we can write 



with 



H' = -Tt (r'(P^)free,free + i^'Wfree.free) 



r 

K' 






^in,out-D 




? 



(3.65) 

(3.66a) 
(3.66b) 



3.2.6 Reduced constraints 

Now we rewrite our constraints, Eqs. (3.9), in terms of the w coordinates. We will 
keep only the parts of the constraint equations that are not automatically satisfied 
because of the ground-state modes. For the expectation value constraints, from Eqs. 
(3.31) and (3.51) we have 



Qout, free (W)fi.ce,freeQout, free ~^ 2^""^'^^^°^^'SS 



out,out X'^'/out.out 

(P) 

out, out (■^) out, out ^free,out {^w} fTee,fTee^ fYee,out 

= £''^(P^)mid,mid^' + const 



const 
T . 



Qiiee.out 2 ^gs,out^gs,out 



(3.67a) 
(3.67b) 



where Wmid means the outside elements of Wfree, i-e. wjVi^+i . . . W2Nin- 

These expectation value matrices must be the same in the vacuum as in the 
vacuum-bounded state. We can accomplish this by requiring that (P«,)mid,mid and 
(W)mid,mid ^rc the same as in the vacuum. Thus we have reduced the problem to one 
that has only A^in(-/Vin + 1) expectation value constraints, regardless of the value of 

A^out- 

For the energy constraint, Eq. (3.9a), we are concerned only with the renormalized 
energy Tr pH — {0\H\0). Thus the constant term in Eq. (3.53) does not contribute, 
and we can use H from Eq. (3.56). Once again there is no dependence on A^out- 
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3.2.7 Derivation based on inside functions 

We would now like to redo the proceeding calculation in a way which does not depend 
on the number of outside oscillators. Then we can remove the box from our system 
by taking L — > oo and — > oo. It appears that we have used the matrices D and D' 
which have an index that runs from 1 to A^out- However, we have used them only in 
particular combinations. The quantities which we need in our calculation are 

1- -f^in,out-D 

2. D'D 

3. D'D''^^ 

4. D^K,,,,,,,D. 

Each of these quantities is an A^in x A^in matrix, so it is reasonable to imagine that 
they do not depend on A"out in the A^ oo limit. 
We proceed as follows: From Eqs. (3.41) we have 

D = (X)o,t,i„ (3.68a) 

D' = (P)in,out . (3.68b) 

Any given element of (X) and (P) has a smooth limit when A'out is taken to infinity. 
It is just a particular expectation value of a half-line of coupled oscillators, which is 
a well-defined problem. We can express the above items in terms of such elements as 
follows: 

1- ^in,out^ depends only on the first row of D which is (X)jVin+i,in so it is already 
well-defined in the limit. 

2. Prom Eq. (3.16b) we have 

(P)(X) = h (3.69) 
so D'D = (P)in,out(X)out,in = " (P)in,in(X)i„,i„, which does not depend on 

3. We expand (P)(P) = ^Z'^^ Z'^Z'^^ Z'K We insert / = ZQqZ'^ here to get 

(P) (P) = ^Z-^^QoZ-^ = . (3.70) 

Then D'D'^ = (P)in,out(P)out,in = ii^in,in- (P)in,in(P)in,in which doCS UOt depend 

on A^out- 

4. Using Eqs. (3.69) and (3.70) we can write 

{X)K = K{X) = (P) and {X)K{X) = (3.71) 
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so 



— / — (X)in,out-^out,out(^)out,in + (X)in,in-^in,out (^)out,in 

+ (^)in,out-^out,in(^)in,in + (^)in,in-^in,m(X)in,in 
— Kout,outD + (X)iii_in((P)in,in ~ -^in,in(X)in,in) 

+ ((P)in,in - (X)in,in^in,in) (X)in,in + (X)in,ini^in,in(X)in,i„ (3.72) 



and so 



D^K,^t,ontD = |/-(X)i,,i,(P)i,,i,-(P)i,,i,(X)i,,i,+ (X)i,,i,Xi,,i,(X)i,,i, (3.73) 

which does not depend on A^out- 

Thus we can now take N ^ oo with A^in fixed and have a well-defined problem in 
terms of K with a finite number of free parameters. 

3.2.8 Calculation of the reduced vacuum 

We are trying to compute K in Eq. (3.54b) in the A^out hmit. We will keep 

iVin and the oscillator spacing Li = L/{N + 1) fixed. With the normalization in Eq. 
(3.12) the vacuum normal mode matrix is given by 



with 



and 



Thus 



I 2 sinfe^/i 



2(A^ + 1) . K 2 . K 
u;^ = y = — sm — . (3.76) 



a=l 



Now we use 



to write 



Li ^ sin kgfi sin kgU 
2(7V + l) A. sin(A;j2) ■ ^ ■ ^ 



cos(^ - 0) - cos(^ + 0) = 2 sin ^ sin (3.78) 



_ Li ^ cos kai/i - i/) - cos A;a(// + ly) , . 
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Using Eq. (3.78) again, for any number a we can write 



cos ka{a - 1) - cos kaa „ . , . 

— = 2sinA;a(a — 1/2) 



sm{ka/2) 



(3.80) 



and thus 



cosM/^- ;-) -cosM/^ + ^) ^2 V sin Ma -1/2). (3.81) 

Sin(/c«/2) a=n-u+l 

If we put Eq. (3.81) into Eq. (3.79) and bring the sum over a inside the sum over a 
we get a sum that we can do, 



. 7ra(a — l/2) 7r(2a — 1) , 

2 E , / = cot , + (-1)" . 



a=l 



N+1 



4(7V + 1) 



(3.82) 



Now we sum this over a. Since /i — u and /i + i/ have the same parity, the (—1)" term 
does not contribute and we get 



Li 



cot 



7r(2a- 1) 



4(iV+l) 4(iV+l) 

^ ' a=iJ,—i/+l ^ ' 



(3.83) 



The sum over N has been eliminated. In the N ^ oo limit, the argument of cot 
goes to zero and so we can use cotx — l/x + 0{x~^) to get 



a=ix—i>+l 

The sum can be done using special functions: 



ijj i IJ, + u + -] - ^|J^|J, — u+- 



where ip is the digamma function iIj{x) = r'{x)/r{x). 

To compute (P) we write the inverse of the normal mode matrix. 



z-} 



2uj„ 



I / — sin kni' 
\ N + 1 ^ 



and 



(3.84) 



(3.85) 



(3.86) 



1 1 ■'^ 

W/ii^ = 2 (^~^^^~^) ^ = jj— a smkafJ, sin kaiy 

^'^ a=l 



N 



— ; -> sm — Sm/CaUSm/Cal^ . 

^ ' a=l 



(3.87) 
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In this case the sum can be done directly. Using Eq. (3.78) we can write 

= C/x-u — Cn+i, (3.88) 

where 



1 A k 



1 



ALi{N + 1) 
Again we take N ^ Xto get 



7r(2A + l) 7r(2A-l) , , 



(3.89) 



1/1 1 
ttZT V2A + 1 " 2A^ 1 J 7rLi(4A2 



and 

^^^^'^ " (4(/. + j.)2-l ~ 4(/.-^)2-l) ■ ^^-^^^ 

Equations (3.85) and (3.91) give (X) and (P) in the A^out hmit. Using these 

values in the procedure of sections 3.2.4 through 3.2.7 we can compute the matrix K 
numerically for a system with inside length Lin but no outside box. 

3.3 Calculation of the vacuum-bounded state 

Once we have computed K we can go on to look for a vacuum-bounded state. We 
arc working entirely with the reduced coordinates. With 

p (X e-'^^' (3.92) 

and H' as in Eqs. (3.65-3.66) we have one number, /3, and two symmetric A^in x 
matrices, -f^mid,mid ^4id,mid' ^bat we can adjust. The constraints (3.9) involve 
one scalar constraint, for total energy, and two Ni^ x A^in symmetric matrices of 
constraints, for tVmWn and Pwm,Pwn- There are equal numbers of equations to satisfy 
and free parameters to adjust, and so, if we are lucky, we will be able to find a 
solution. If we do find a solution, we know it is unique from the arguments of section 
2.2.3. 

3.3.1 Expectation values 

To actually solve these equations we will need to compute the expectation values of 
WmWn and PwmPwn given the density matrix (3.92). To do this we find the normal 
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modes of H'. The Hamiltonian H' gives rise to classical equations of motion 

^^^-T'X'Wfree, (3.93) 

so we look for eigenvectors wggg that satisfy 

T'K'wg^^ = a;>f^,, . (3.94) 

The eigenvectors will be complete, so that we can define new coordinates Ua via 
Wfree = S'^^aWggg, which wiU then obey the equations of motion 



^-ujiua. (3.95) 



We can choose the norms of the eigenvectors so that^ '^free' -^'^hee ~ ^a^af) and group 
the eigenvectors into a matrix U via Una — w'^. Then we will find that = '^^'^vi 
or 

K' = U-^^n'^U-^ and T'^UU^ (3.96) 
where JIq,^ = ^^a^a/s- We can then substitute 



Wfree = UU (P^)free = ^"'^P« (3.97) 



into H' to get 



= ^E(^- (3-98) 

a a 

This is the Hamiltonian for a set of (fictitious) uncoupled oscillators with frequencies 
cUa- The expectation values of the u and are easily computed, 

= -^5a/3Coth-^ (3.99a) 

{Pu^Pu,) = Y^a^ coth ^ (3.99b) 

SO that in terms of the w coordinates we have 

{WmWn) = TT—^^ocUna COth (3.100a) 

a 



^This is a different normalization than we used for Z in section 3.2.1. In the present normahzation 
the original vacuum modes would all have the same amplitude. In section 3.2.1 the modes were 

— 1/2 

multiplied by a factor of tOa relative to the convention here. 



48 



Chapter 3. The Problem 



{Pn,r.PwJ = E Y^a^^^n COth . (3.100b) 

a 

In section 4.3 we will compute these expectation values numerically in the vacuum 
state, given hy H' = H and (3 = oo, and require that they have the same values in 
the vacuum-bounded state with finite (3. 



3.4 The nature of the results 

Before looking at the results of our computations, we would like to learn as much as 
possible about the form that our answers must take. We will take the frequencies cua 
as given here and look at the form of the normal modes w". 
The vectors w^Jigg satisfy the equation 



(3.101) 



with 



and 



K' 











r 



-^in,m 


-^in,out-C^ 




K' 

mid,mid 



(3.102) 



(3.103) 



However, we have tridiagonalized K, so only the lower left element of Ki^^o^tD is 
nonzero. In fact, in our numerical work we will find that this element is always —g, 
just as it was in the original i^^in,out- Thus we have 



K' 



( K- 
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(3.104) 



Writing out Eq. (3.101) in components, we get 



(3.105) 



Taking only the inside components of the eigenvalue equation, we see that 



K' = u'^w°' 

an n a a 



for a < iVin. That is to say. 



(3.106) 



2gw^ - gw^ = 
-gw'^ + 2gw^ - gw^ = 



(3.107a) 
(3.107b) 



3.4. The nature of the results 
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-9w%.^-i + '^9w%.^ - 9w%.^+i = i^a^N,^ ■ (3.107c) 

Taking Ua fixed, there are A^in equations involving Ni^ + 1 unknown components of 
w". However, the equations are invariant under a uniform rescahng of w". Thus these 
A'in equations fix for = 1 . . . A^in + 1, except for normalization. The equations 
are readily solved, and the solution is 

w^^Ksmfik'^, (3.108) 

where 

cosA;; = l--^, (3.109) 

and N'j, is an unknown normalization factor. Here k'^ and A*"^ can be complex, but 
w" must be real. When k'^ is real we will write ka = k'^ and Na = N'^ and call this a 
"normal" mode. When A;^ is complex we can write 

k'^ = 7r + ika (3.110a) 
N'^ = iN„ (3.110b) 
= {-y-^N^ sinh ^Jika (3.110c) 

with ka and Na real. We will refer to these as "abnormal" modes. 

A similar calculation can be done for U~^, the inverse of the eigenvalue matrix U. 
In this case we will find that 

U-l=Ksmak'^. (3.111) 

Equation (3.111) has the same form as Eq. (3.108), but applies only for a = 1 . . . N[^. 
That makes Eq. (3.111) less useful than Eq. (3.108) for establishing a connection 
between the inside and the outside region, and we will not use it further. 



Chapter 4 
Numerical Solution 



Two numerical calculations are necessary to solve the problem. First we must fol- 
low the procedure of sections 3.2.4-3.2.8 to find the reduced Hamiltonian H and its 
vacuum state. Then we must search for parameter values T^^^j -^mid mid ^^"^ 
which produce a state with a given energy Eq but the same expectation values of the 
outside oscillators. 

In fact, we follow a slight variant of the above plan. Instead of fixing Eq, we 
hold (3 fixed. The resulting state has some Eq and solves the problem of maximizing 
the entropy for that £"0, whatever it is. By varying /5 we can find states for various 
energies. 

In the end we do not use these results directly to calculate bounds on the entropy. 
Instead we derive a general principle from the numerical calculations, and use this 
principle as an ansatz to derive a bound in the next chapter. 

4.1 Numerical procedures 

Finding the reduced vacuum is a fairly straightforward problem in numerical analysis. 

The number of steps grows as Nf^. However, in order to produce accurate results 
for iVin ^ 6 it is necessary to use very high-precision floating-point numbers. The 
necessary number of bits of mantissa in the representation appears to be about lOA^in. 

After finding the reduced vacuum we need to solve a set of simultaneous nonlinear 
equations. Such problems are in general quite difficult, and require an iterative search 
for the correct parameter values. Here at least we know from section 2.2.3 that there 
cannot be more than one solution. Although we have not been able to prove that a 
solution always exists, in the numerical work we have always been able to find one. 

Once the search is sufficiently close to the correct answer it is possible to use 
Newton's method, which converges quadratically, i.e., the number of correct digits 
doubles every step. However, the basin of attraction for Newton's method can be 
quite small and difficult to find. When Newton's method does not work, it is necessary 
to use some other procedure to make progress toward the solution. Here we used 
the Powell hybrid method [37,38]. This method moves in a direction which is a 
combination of the direction Newton's method would suggest and the direction of 
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steepest descent in the mean square error of the function values. Such a method has 
the difficulty that it can get stuck at a local minimum of the mean square error that 
is not a solution. We have been lucky in that there do not seem to be such local 
minima in this problem.^ 

In the case that there are no local minima of the mean square error in the func- 
tion values, Powell's method is guaranteed to converge from any starting point [37]. 
However, it often converges quite slowly for large systems, requiring many thousands 
of iterations to make progress. This has limited our numerical solutions to problems 
with iVin < 15. 

The code was written in Lisp and executed on DEC^^'' Alpha™ workstations. 
All results presented here were computed using at least 38 decimal digits of precision. 
The most precise calculation (the reduced vacuum for A^in = 175) used 636 digits. 
The parameters found in Powell's method reproduce the desired expectation values 
to at least 17 significant digits. 

We have made use of many routines from Numerical Recipes [2]. However, the 
code that implements Powell's method was written from scratch in Lisp following the 
outlines of [37,38] and has many extra features including dynamic increase of working 
precision as the solution converges. 



4.2 Reduction of the vacuum 

First we follow the procedures of sections 3.2.4-3.2.8 to calculate the matrix K which 
gives the ground state of the Wfree coordinates. The result is, of course, a tridiagonal 
matrix which gives a set of self-couplings and nearest-neighbor couplings for the fic- 
titious oscillators Wfree- We can express these couphngs as multiples of the couphngs 
for a regular chain of A^free oscillators with spacing Li. Thus we write the self-coupling 
as 

K,, = -9U {t,) (4-1) 
and the nearest-neighbor coupling as 

K,,,+i = 2^/^+1/2 {t,) ■ (4-2) 

These coupling coefficients converge rapidly to a universal form f{x) where Li// — > 
X in the continuum limit. Some results are shown in Fig. 4-1. We can see that 
f{x) ~ 1 until a; ~ 1.4 at which point it begins to fall and asymptotically approaches 
as a; — > 2. For values of x near 2, f{x) is well fit by 

f{x)^a{2-xy (4.3) 



^The current formulation of the problem has 0{N?^) parameters and equations. It is possible 
to use the normal-mode frequencies and normalizations as our parameters, which gives a problem 
with only 0(A''in) parameters and equations. However, this reduction loses the property of having 
no local minima, and so in fact makes the problem harder to solve. 
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Nin = 20 

Nin = 50 

Nin = 100 



1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

X 

Figure 4-1: The ratio of the coupUng coefficients in the reduced prob- 
lem to what they would be in a regular problem of Nfree oscillators, 
plotted against x — Li/i, for L,n = 1. 



with 3.2, as shown in Fie;. 4-2. 

Some typical normal modes of the reduced vacuum are shown in Fig. 4-3. They 
are sine waves while f{x) ~ 1 and then begin to oscillate faster and faster as f{x) 
shrinks. At first the amplitude of the oscillations grows but for larger x it shrinks 
rapidly to zero. The wavenumbers in the inside region (and thus the frequencies) are 
smaller than we would find for a rigid box because most of the oscillations are in the 
part of the outside region where f{x)<^l. In fact, as — > cxd we would expect the 
low-lying frequencies to go to zero, for the following reason. 

We can find the frequencies by computing the normal modes of a Hamiltonian 

[ dxdy {T{x - y)7v{x)n{y) + K{x - y)<j){x)<j){y)) , (4.4) 
which requires solving the eigenvector equation 

/ dydzT{x - y)K{y - z)g{z) = Xg{x) (4.5) 
Jo 

with the boundary conditions 



g{0) = 
giL) = . 



(4.6a) 
(4.6b) 
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data <> 
fit — - 



oj>oooooo6ooooo 



1.8 



1.82 1.84 1.86 1.88 1.9 1.92 1.94 1.! 

X 



1.98 2 



Figure 4-2: The coupling cocfRcient ratio for Ni^ — 175 in the region 
X > 1.8 and the fit f{x) = 3.2(2 - x^. 




Figure 4-3: The 5 and 25 normal modes in the reduced vacuum, 
computed with Ni^ — 175. 



4.3. Computation of the vacuum-bounded state 
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Since Eq. (4.5) is a second-order differential equation we expect two degrees of freedom 
in the solution. However, one degree of freedom is manifestly the overall scale, which 
does not affect the boundary conditions. Since there are two boundary conditions but 
only one free parameter, we can expect to find solutions only for particular values 
of A. For example, for the usual scalar field Hamiltonian the general solution to Eq. 
(4.5) would be 

g{x) = csm{\/\x + S) . (4.7) 

To satisfy Eqs. (4.6) we need to choose 6 = and \/A = mr/L for some integer n. 
However, if we use 




with f{x) a{2—x)^, as suggested by Fig. 4-2, we will get a continuum of frequencies. 
The problem is that since /(x) — > as x — > 2 the boundary condition there does not 
really constrain g{x). There can be arbitrary changes in g(x) near x = 2 and so Eq. 
(4.6b) can always be satisfied. Since there is only one effective boundary condition 
and one effective degree of freedom, we expect to be able to find a solution for any A. 
Thus in the continuum limit there are modes with arbitrarily low frequencies. This is 
not an unreasonable conclusion, since although the range of x is finite, we are using 
it to represent the infinite half-line. In the infinite vacuum there is no right-hand 
boundary condition, and there are modes of every frequency. 

This conclusion is confirmed by numerical results. In Fig. 4-4 we plot the lowest 
normal-mode frequency versus A^in. As shown in the figure, the frequencies are well 
fit by a curve 

ax-^ - x-'^ (4.9) 

with a ~ 1.81. If this form is correct, in the A^in — > oo limit the lowest frequency goes 
to zero. 

If we go to a vacuum-bounded state we will introduce some finite temperature. We 
then expect that the non-zero temperature will increase the frequencies in such a way 
that there are only a finite number of low-lying modes and thus a finite entropy. How- 
ever, in the limit where T ^ we do expect an entropy-to-energy relation equivalent 
to a system with infinitesimal frequencies. This is discussed in the appendix. 



4.3 Computation of the vacuum-bounded state 

Once we have computed K and the vacuum expectation values of WmWn and 

we can proceed to look for the vacuum-bounded state at a particular (3. To do this 

we vary the A^in(-/Vin + 1) independent components of T^id^mid -^mid,mid 

those which reproduce the same expectation values of the A^in(-/Vin + 1) independent 

operators WmWn and PwmPw„-i given in Eqs. (3.100). 

To understand the numerical solution we look at the normal mode frequencies 
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Figure 4-4: The lowest frequency of the reduced vacuum and the fit 

and the forms of the normal modes. When the mode is "normal" (i.e. real k'^ in Eq. 
(3.108)) there is a sine wave in the inside region. When the mode is "abnormal" there 
is essentially a growing exponential. Typical modes for a small number of oscillators 
are shown in Figs. 4-5 and 4-6. 

As N becomes large, each "normal" mode and its frequency smoothly approach 
a limit, providing that we use a normalization appropriate for the continuum, which 
means that each mode must be rescaled by . In Fig. 4-7 wc show the first normal 
mode for various values of A^in. Note that this mode docs not appear to come down 
to zero ai X = 2. As discussed in section 4.2, this happens because f{x) is going to 
at X — 2 and so there is not really any coupling to the boundary at that point. 

As A^ increases, each abnormal mode and its frequency undergo a smooth evolu- 
tion, until at some point it disappears from the set of abnormal modes and is replaced 
by a normal mode with very similar form in the outside region, as shown in Fig. 4-8. 
Because of this behavior, we believe that if we could solve the continuum behavior 
directly we would find just the "normal" modes. 



4.3.1 Evenly spaced wavenumbers 

To address the problem directly for large energies would require numerical solutions 
for large numbers of oscillators, which is computationally intractable. Instead we 
would like to extract from the computations in the accessible regime a statement 
which will allow us to extend our arguments to larger energies. The most striking 
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Figure 4-5: Modes and frequencies for the "normal" modes of a 
system with Li^ — 1.0, P — 2, Ni^ — 3, A^free — 6. 
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Figure 4-6: Modes and frequencies for the "abnormal" modes of a 
system with Li^ — 1.0, P — 2, Njn — 3, A^free — 6. 
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Figure 4-7: The first normal mode, rescaled to the continuum nor- 
malization, for various different numbers of oscillators, with Lin = 
1.0, P = 0.5. 



such result is that the wavenumbers of the "normal" modes are nearly evenly spaced. 
The larger the energy, the more accurate is this approximation. In Fig. 4-9, we 
show the wavenumbers for Lin = 1-0 and f3 = 0.5, which give an energy of about 
0.094. Even at this low energy the fit is good to within a few percent of the typical 
wavenumber. For larger energies the points will lie correspondingly closer to the line. 

Since we have set N and L to oo, our problem has only two dimensionful pa- 
rameters, Lin and Eq. Thus there is only one dimensionless parameter, L^-^Eq, that 
characterizes the problem. In the 3-dimensional black hole problem the equivalent 
parameter is REq, which for the parameters in Eq. (1.3) is about 10^^. Thus for 
application to black holes we are interested only in very large values of Li^Eo, for 
which the linear approximation for the "normal" modes will be very good. 

Since the "normal" mode wavenumbers extend up to A; ~ tt, the number of "nor- 
mal" modes will be given an integer A^norm ~ where ki is the interval between 
wavenumbers. In addition there will be iVfree — -/Vnorm "abnormal" modes, with fre- 
quencies cUa > 2/Li. For N^j^ ^ Lin//3 these modes do not contribute to the entropy, 
because they are exponentially suppressed. 
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Figure 4-8: An abnormal mode for Li^ = 1.0, (3 — 0.5 and various 
values of Nj^. The normalizations are chosen to be similar in the 
outside region. As increases the mode changes smoothly until 
at Njn = 12 there is no corresponding abnormal mode, but instead 
one of the new normal modes has a very similar form (but different 
normalization ). 
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Figure 4-9: The numerically computed wavcnumbcrs compared with 
the best-Gt line through the origin for Lin — 1-0, /3 — 0.5, Njn — 12. 



Chapter 5 

The Entropy Bound 



As discussed in section 4.3.1, wc arc interested in vacuum-bounded states for quite 
large energies. When Eq is large many different modes contribute to the entropy. 
To get an accurate result for a system in which many modes are important requires 
using many oscillators, which is computationally intractable. Instead, we take the 
even wavcnumber spacing of section 4.3.1 as given and derive a bound on the entropy 
from that ansatz. 

5.1 The first outside oscillator 

With the wavenumber spacing fixed we have one free parameter, the spacing ki, 
which depends on the energy Eq. To fix ki we we examine the expectation value of 
wj^. Since this is an outside operator, it must have the same value in the vacuum 
as in a vacuum-bounded state. The vacuum value is straightforward to compute, 
and we show below that in the vacuum-bounded state the value depends only on the 
wavenumbers. In fact, we will be able to derive only an upper bound in the vacuum 
case and only a lower bound in the case of a vacuum-bounded state. However, these 
bounds are sufficient to derive a lower bound on ki and thus an upper bound on S{E). 

5.1.1 The vacuum 

In section 3.2.8 we computed the values of (X). To convert to w coordinates we 
proceed as follows: Prom section 3.2.4, x = Qw and so 



out 



Dw^id + >^out,gsW, 



gs • 




In the numerical work we found that with our choice for D we got 




(5.2) 
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Since 



in, mid 



-^in,out-D 







-9\ 



D 



(5.3) 



it follows that the first row of D is (1, ... 0). Thus from Eq. (5.1) we get 

Xn,^+1 = WiVi„+l + (5^out,gs)l • Wgs (5.4) 

where (yout,gs)i denotes the first row of l^ut,gs- Now Wgs,gs = (1/2)/ and Wgs,free = 0, 
so 



(^ATin+l) — ('"^ATin+l) + 2 (■^gs,out^gT,out)ii 

Since the last term is non-negative, we have 
Prom Eq. (3.85) we have 



2 \vac 

■-Afin + l/ ' 



+1/ 



hi 

27r 



(5.5) 



(5.6) 



(5.7) 



We would like to evaluate this expression in the Ni^ — > oo limit with fixed. There 
is a pref actor of Li, which goes to zero in this limit, but that is just an artifact of the 
conventions we have used for the discrete problem, and will appear in the finite-energy 
vacuum-bounded states as well. For large x. 



il^^x) ~ Ina; -|- 0(l/a;) , 



(5.8) 



so without the prefactor there is a logarithmic divergence. We are interested in the 
In A^in term, and in the constant term, but we will ignore any terms of order 1/A^in or 
lower. 

We use Eq. (5.8) and %l}{l/2) = —7 — 2 In 2, where 7 is Euler's constant, to get 



2n 



In 2Nin + 'y + 2hi2 + 



and so 



27r 



ln8^"in + 7 + 



(5.9) 



(5.10) 



5.1.2 Evenly spaced wavenumbers 

Now we will compute the same correlator in the vacuum-bounded system, using the 
ansatz that the wavenumbers are multiples of some spacing ki. Prom Eq. (3.100a) we 
have 

(5.11) 
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There are A^norm ~ 7r//ci normal modes which are sine waves in the inside region. 
These modes are equivalent to the modes that we would have for a problem with a 
rigid boundary at position 

4n = A^norm^l ~ ^^1 ■ (5.12) 
Ki 

There are also A^free — -^norm abnormal modes. We do not know how to compute the 
contribution of the the abnormal modes to the correlator. However, the contribution 
from each mode is positive, so by taking only normal modes we will find a lower 
bound on {w%.^^i). 

Prom Eq. (3.108) for the normal modes we have 

w'^ — NnSin jikn (5.13) 

for = 1 . . . A^in + 1 . Putting this in Eq. (5.11) we get 

(w^ +i) = / coth I —ujn ) sin^ kn{Ni^ + 1) + abnormal modes 

^norm f Q ^\ 

J] ^ (f + 1) (^-1^) 



> 



since the abnormal mode contribution is positive. 

The important point here is that wc know w'^ for /x up to A^i^ + 1, and we know 
{w^w^) for n and v down to A^in + 1. Thus taking /i = u = A^in + 1 gives the unique 
correlator for which we know the components that go into the expression for the 
correlator while also knowing that the correlator must have the same value as in the 
vacuum. The same argument does not work for (P^^. +iPwn. because we know 
the inverse mode matrix only up to /i = A^in and not fi = A^in + 1. 

Now we will compute the right-hand side of Eq. (5.14) in the limit where A^in and 
A^norm are large. We will ignore all terms of 0(1/A^in) or 0(l/A'norm) and so take 

-^'in/A^norm = AnMn- We USe 

TTTi 

kn = nki = — (5.15) 

'norm 

and 

2 . kn 2 . Trn 
= — sm— = — sm— . (5.16) 

Pirst we must derive the normalizations of the normal modes. With our choice of 
normalization. 
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Thus if a is an inside oscillator, then for any /x, 

Y.«-^^^ (5.18) 

n 

so 

-^norm 

<^a/i — -^n ^nfl sin + abnormal modes . (5.19) 

n = l 

We would like to use Eq. (5.19) to determine the normalization factors Nn, but first 
we have to dispose of the abnormal mode term. Later we will find that A^norm is close 
to iVin. Using this, and since we are working in the limit where iVin is large, we can 
choose a such that 

K a < 2Ari^ - Nr,orm ~ A^in . (5.20) 
Each abnormal mode n contributes 

sinh kna sinh A;„/x = 6^^^ (5-21) 

to the right-hand side of Eq. (5.19). However, this same mode contributes 

Li 2 sinh^fcn(iVm + 1) 

4 " cosh(A;„/2) ^ ^ 

to the sum for in Eq. (5.14). Since (fw^j^+i) < Li/(27r)[lnA^in + 0(1)] it 

follows that 

^2sinh2/c„(iVi„ + l) ^ 2 



C0sh(fc„/2) TT 



for each n and thus that 



(„) 2 sinh /c^a sinh A;„/i cosh(A;„/2) 

<^abn < ■ , 2, ... ^.x In . (5.24) 

TTSmh Kn[-rJ in + 1) 

The exact values of the A;„ for the abnormal modes vary with A^in- For successive 
values, each A;„ decreases toward until the corresponding mode converts to 
a normal mode as described in section 4.3. By avoiding the points where these 
"conversions" are about to take place it is possible to find a sequence of values for 
A^in which have all A;„ ^ 1. For such values, 

< g-fc„(2iVi„-o-(U+3/2) _ (5.25) 

From Eq. (5.20) the exponent in Eq. (5.25) is <^ —1, so S''^^ is exponentially small. 
Since there are at most A^in abnormal modes, their total contribution to (5.19) is 
exponentially suppressed. 

Thus we ignore the abnormal modes in Eq. (5.19), multiply by sin/cmA* and sum 
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over jj, to get 

-^norm -^norm 

sin kmd = sin kmlJ' sin A;„a sin /c^/x 

//=1 n=l 

orm rm 

= E ^nSinA;„a E sin ™ sin ^ 

n=l ;U=1 
orm 7\ 7- 

sm A;„a ^ — ^rnn 

n=l 

= <%^sinfc^a, (5.26) 
from which we conclude that 

< = T7^ (5.27) 

-' 'norm 

for each mode m. Putting this in Eq. (5.14) gives 

{wn,^+i) ^ 2^ iTr coth [-uJn] . (5.28) 

Note that Eq. (5.28) does not depend on but only iVnorm- 

Let us define a dimensionless parameterization of the departure from the vacuum 
state, 

TV — TV- L'- — L- 

-^norm -^in 

r' ^ UiJ'^U-jp. (5.29b) 
Then sin^ /Cn(^in + 1) = sin^ n7r(l — A) = sin^ nvrA. With uj^ from Eq. (5.16) we get 

orm • *^ \ / O \ 

^ Y E AT sin^ ^"^^ 2A^) • ^'-'^^ 

In the limit of large A^frcc- the argument of coth becomes /37rn/(2LiA^norm) ~ 
7rn/(2r'). We expand cothx = 1 + 2/(e^' — 1) to get 

2 \ Lx^^^ sin^riTrA / 2 \ /^o-,\ 

We work first with the term not involving r'. We expand the numerator using sin^ x ~ 
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(1 — cos2a;)/2 to get 

^ ( ^ ""-""^ 1 . (5.32) 

4 I TV sin N sin / ^ ^ 

^ 1 \ ' norm "^^^ 9 at ' norm ''^J^i 9 / 

This first term is the one that diverges as A^norm —>■ oo. We can separate out the 
divergent part to get 

orm / 1 O \ 

h y 1 + ^ y 1 A . (5.33) 

2^ ^ ^ 4 ^ l^normSin^ nU ) ^ ^ 

The first term of Eq. (5.33) gives 

^(ln7Vnorm + 7) ■ (5-34) 
ZTT 

The second term of Eq. (5.33) is finite and can be converted to an integral in the 
-/Vnorm OO fimit, to give 



4 Jo Vsin^ TTX. 



The remaining term of Eq. (5.32) is 





In 




nx\ 




tan — - 


1 27r 




\x 


4 J. 



^Mn - . (5.35) 



27r TT 



iWn cos2n7rA 
-T^TV sin^H^ - ^^-^^^ 

^ „ 1 ' norm ^^J^^ o^^ 

To compute this we use 1/ svax — cscx — 1/x + x/6 + 7x^/360 + • • • to get 

The first term can be summed in the A^norm oo fimit, 

1 ■^cos2n7rA 1 , . * x 

-— > = — ln 2sin7rA . 5.38 

27r^ n 27r ^ ^ ^ ^ 

The rest of the terms do not contribute. Because of the osciUations of the cosine, 
Xli ™ cosnTrA goes as N^orm rather than N^orla thus is killed by the corre- 
sponding iV^+^ in the denominator. 

Putting Eqs. (5.34), (5.35) and (5.38) together, the first term on the right of Eq. 
(5.31) gives 

Ll f SiVnorm sin TT A \ 

-fin + 7). (5.39) 



5.1. The first outside oscillator 



67 



We now look at the second term of Eq. (5.31), 



siir iiTvA 1 

n=l ^jVnorm ^ ' 



Here the divergence is cut off by the exponential in the denominator. We again 
expand using l/sina: = l/a; + a;/6H — •. The first term has no iVnorm dependence and 
we can extend the sum to oo. In the next term, the sum is cut off by the exponential 
in the denominator, leading to a term of order (r'/A^norm)^- Further terms have higher 
powers of r'/A^norm- In the limit A^norm — > oo we ignore all these terms, which leaves 



2Li sin"' rnrA 



TT ie^^h' - 1) ■ 



(5.41) 



We are interested in the high-energy limit, for which r' ^ 1. Later we will see that 
A is of order In r'/r' -C 1. Thus the summand in Eq. (5.41) is slowly varying and we 
can convert the sum into an integral, 

2Li /""^ sin^ TT Ax (5 42) 



TT Jo X (e'^^/^' - 1) ■ 

The error in Eq. (5.42) is approximately the term that we would have for n — > in 
Eq. (5.41). Taking this limit we find that the error has order A^r' ~ (In r')^/r' ^ 1, 
so our approximation is good. The integral in Eq. (5.42) can be done and the result 
is 

Lir'A+^lni-^— ^. (5.43) 



27T inr'A 
Putting together Eqs. (5.39) and (5.43) we find 

/ 2 \>^lA 8iVno.n,sin7rA , l^-4-r^\ 

(^iv.+i) ^ ^ (^In + 7 + 27rr A + In j ■ (5.44) 

Now we set (w^.^^i) = from Eq. (5.10) to get 

87V- 87V sin ttA 1 — p-^^r'A 

In ^ > In ^^^normSmTTA ^ ^^^/^ + lu ^ , (5.45) 

or 

2nT'A + In ^—^—^ < In ^ . (5.46) 

47rr'A TVnormSmTrA ^ ^ 

Now we use TVin/TVnorm = 1 — A from Eq. (5.29a) and approximate sIuttA ^ ttA since 
A is small. Since this is already 0(A) we then approximate (1 — A)/Apa 1/Ato get 

1 — e""^'^'^'^ 1 
27rT'A + In ^ ^. < In — . (5.47) 
47rT'A ttA 
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Thus A < Amax where 

^^l^T^^^T^- (5.48) 

If instead of r' = L\J^ we use 

T = LinT (5.49) 

we will make an error of order Amaxj which wc expect to be small. We ignore this 
second order contribution and take r' as r in Eq. (5.48), 

1 4r 

(5.50) 



""^ 27rT 1 - 
If we ignore e"^'^'^^™'^ in the denominator we get 

1 



zvrr 



ln4T. (5.51) 



Using this we find that e ^^'^^'^^^ = (4r) ^ ^ 1 since r ^ 1, which justifies ignoring 
this term. We will also ignore In 4 by comparison with In r. Thus we conclude 

and 

4n < i^in + ^ In Li„T + C . (5.53) 

The equivalent system is larger by at most a thermal wavelength times a logarithmic 
factor depending on the inside size. 



5.2 Propagation of bounds 

In the previous section we derived an expression that gives the frequencies, and thus 
the entropy, for a vacuum-bounded system at a given temperature T = 1/(3. Given 
such an expression, we would like to compute the entropy as a function of energy. 
Unfortunately the energy is not simple to compute from the frequencies alone. ^ How- 
ever, we can easily compare the entropy of the vacuum-bounded system to that of 
a system with a rigid boundary at Lin and the same temperature. To make this 
comparison at fixed energy instead, we proceed as follows. 



^Such a computation can be done, but since E needs to be renormalized against the ground-state 
energy of the entire system, the result depends sensitively on the frequencies and normalizations 
even for very high-energy modes. 
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Consider the free energy F = E — TS which has dF = —SdT. Integrating gives 

E-TS^- [ S{T')dT' . (5.54) 



^0 

Let 5"^^ and T^^ denote the entropy and temperature of the system with a rigid 
boundary at Lin, and S and T denote those for the vacuum bounded system. For any 
quantity A let 5A{T) denote the difference between vacuum-bounded and rigid box 
systems at fixed temperature, 6A{T) = A{T) — A^^{T). and 6A{E) denote the same 
difference at fixed energy, 5A[E) = A{E) — A'^^^E). With E fixed we compare the 
differences (to first order) in the two sides of Eq. (5.54) between the vacuum-bounded 
and rigid box systems, 

-TSS{E) - ST{E)S = -ST{E)S - I SS{T')dT' , (5.55) 

Jo 

where the first term on the right-hand side comes from the change in the integration 
hmit. Thus 

1 

SS{E) = - 5S{T')dT' . (5.56) 
T Jo 

5.3 The final entropy bound 

Now we apply Eq. (5.56) to the case of section 5.1 where 



n 



— and L'in- Li„(l + A) (5.57) 



-^in 



with 

ALin < -i-lnLinT. (5.58) 
zttJ 

At any given temperature, the vacuum-bounded system has the entropy iS'(T) of 
a system of length L'n. Now in a one-dimensional system the entropy density is 
proportional to the temperature, 

,S"-^ = |Linr, (5.59) 
and thus the entropy difference between vacuum-bounded and rigid box systems is 

bS{T) = I ALi,T < i In UJ" . (5.60) 
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Using Eq. (5.56) we get^ 

SS{E) < 1 i In L,XdT' = ^ (In L,^T - 1) . (5.61) 

Since we are ignoring terms of order 1 by comparison with those of order In T, we can 
write 

SS{E) < ^ In LinT f« ^ In S'"" . (5.62) 

Thus we conclude that the vacuum-bounded condition closely approximates the 
rigid box of length L^n- For the same energy, the vacuum-bounded condition allows 
slightly more entropy. The entropy difference grows at most logarithmically with rigid 
box entropy. For high energies, S"^^ ^ 1 so we conclude that 6S < S'^. 



5.4 Discussion 

We have introduced a new way of specifying that matter and energy are confined 
to a particular region of space. Rather than giving a boundary condition per se, we 
specify a condition on a density matrix describing the state of the overall system. 
We require that any measurement which does not look into the inside region cannot 
distinguish our system from the vacuum. This avoids certain difficulties such as the 
Casimir energy that results from the introduction of a boundary and the geometric 
entropy [31,32] that results from ignoring part of a system. For these vacuum-bounded 
states, we consider the problem of finding the maximum-entropy state for a given total 
energy. This is analogous to the problem of finding the thermal state in a system 
with a rigid boundary. 

Unfortunately, the vacuum-bounded problem is more difficult than the analogous 
problem with a rigid boundary and we must resort to working in one dimension and 
to numerical solution on a lattice. It is, however, possible to reduce the problem to a 
finite number of degrees of freedom, even when the outside region is infinitely large. 
From the numerical solution we justify the ansatz that the continuum wavenumbers 
are evenly spaced in this problem. Using this ansatz we compute an upper bound on 
the entropy of a vacuum-bounded state, and show that for high energies {ER » 1) 
the entropy approaches that of a system with rigid boundaries. Of course this is what 
one would expect for a system whose typical wavelengths are much shorter than the 
size of the inside region. 

To apply this result to an evaporating black hole we look at the state produced 
by the black hole after evaporation [5] . Since our calculation was one-dimensional we 
must assume that the similarity between the vacuum-bounded state and the thermal 
state with a rigid boundary extends to three dimensions. Then we infer that very 



^It happens that SS{E) and 6S{T) are approximately the same, but that is a particular property of 
the system at hand. For example, if A were a constant we would have SS{E) = 1/T cALi^T'dT' = 
2cALi„T = 2SS{T). 
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little entropy can be emitted in the final explosion, in accord with the results of 
Aharonov, Casher and Nussinov [24] and Preskill [5]. For example, a black hole 
formed in the big bang with mass of order lO^^^i would be evaporating today. During 
its hfe it would have radiated entropy S ~ 10^^. Now wc assume that the entropy 
of the final explosion has energy E ~ lO^^erg contained in radius R ~ 10~^^cm as 
in section 1.2.1, and that the maximum entropy is not too different from that of 
a spherical box, in accord with our one-dimensional result. Then we find that the 
final explosion can emit only entropy S ~ 10^°, which is a factor of 10^^ less than 
what was emitted earlier in the thermal radiation. The choice of Tunk is somewhat 
arbitrary, but whatever value one chooses there is some fixed bound on the emission 
of entropy after T-^nh is reached. By considering a sufficiently large starting black 
hole, and thus sufficient entropy emission at early times, one always finds that the 
late time information is too little to produce a final pure state. 

This argument means that a black hole must not evaporate completely but rather 
leave a remnant or remnants, that information must be lost, or else that the Hawking 
radiation is not exactly thermal, even at very early times [16]. 
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Appendix 

Low Energy Results 



In the low-energy regime we do not see the hnear wavenumber relation that we 
see in the high-energy case (Fig. 4-9). Instead, for sufficiently low temperature, the 
wavenumbers and frequencies are nearly the same as in the reduced vacuum. These 
frequencies are much lower than in a rigid box with the same Lin. At low temperatures, 
the entropy depends only on the low-lying frequencies and on j3. Thus we expect that 
there will be significantly more entropy in a vacuum-bounded state than a rigid box 
state of the same energy. While we don't know how to construct an analytic proof of 
this claim, we will outline a general argument here, and make a conjecture supported 
by numerical data. 

In the low-energy regime we can make a first-order expansion around the vac- 
uum. To do this we note that the only dependence on f3 in our equations is through 
coth(/3a;a/2) in Eqs. (3.100). For large /? we can approximate 

coth ^ 1 + 2e-^'^« . (A.l) 

The change in coth{(3u!a/2) is the largest for the smallest frequency, which we will 
call uji. We will ignore e~^'^" for larger uja by comparison with e~^'^^. Thus we take 

5coth^ = 2e-^'^™^ = 2e (A.2a) 
(5coth ^ = for a > 1 . (A.2b) 



Then we write 

-^mid,mid ~ ^mid,mid ~l~ <^^mid,mid (A. 3a) 

-^mid,mid ~ ^mid,mid + <^-^mid,mid (A. 3b) 

where 5i^mid,mid and 5Tmid,mid are 0(e). These changes give rise to 0(e) changes in U 
and the cUa, which in turn give rise to 0(e) changes in {wmWn) and {PwmPw„)- 
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Since overall {wmWn) and {PwmPw„) cannot change we must have 

= S{WmWn) = ^ 
a 



= S{P^^P^J = ^ ( 

a 

We thus have iVin(A^in + 1) linear equations for A^in(A^in + 1) unknown values of 
5Tjnid,mid and (5-ft'mid,mid5 which are readily solved. Since the inhomogeneous part of 
these equations is 0(e), all the results must be 0(e) as well. In particular, the 6uJa 
are 0(e). Now if T is very small as compared to all the a;^, then e will be small as 
compared to all the parameters of the problem, and so the first-order approximation 
will be good. For any fixed number of oscillators A^in there will be some minimum 
frequency uji, and if we take /5 <^ l/cui we will always be in this regime. 

Now the entropy S depends only on /3 and the cUa- Since the modes are uncoupled, 

>S = J]>Si(/3u;«) (A.5) 

a 

with 

S^{Pu;) = - ln(l - e"^-) + . (A.6) 

Since (3u!a 1 all the terms are very small, and the cui term dominates, 

S ^ S,{u;,) = (1 + /?wi)e-'^"i + 0(e-2'^"i) . (A.7) 

Since e drops exponentially with increasing /3, we expect that for (3 large enough, 
PScui ^ 1 so that 

S=il + PuDe-^'^^^^ + 0(e2) . (A.8) 
The value of S given in Eq. (A.8) is the one we would get from a rigid box with length 

4n = 7r/<'- (A.9) 

To approximate the energy, we proceed along the lines of section 5.2. The direct 
calculation is made difficult by the fact that, while H' differs from H only by 0(e), 
we must subtract from both Hamiltonians a large ground-state energy. Instead we 
work by integrating on T. Prom Eq. (5.54) we have 

E{T) = TS{T)- [ S{T')dT' . (A.IO) 

Now L'in depends only on u'^, which depends on iVin but not on /3. If Eq. (A.8) is 
valid for a particular N-,^ at /3 = 1/T is it valid for T <T and (3' = l/T' > /3. Thus 



'^^a Tyvacry-vac J \ grj Tjyac : ]_7y !tt \ 



+-4t^mTt^nTe (A.4a) 

TT^^'^~ rrvac" i a sjj—ljj\ac~ i a Tjva,c~ K[J~^\ 

2 am an ' 2 am an ' 2 am an J 

+u;rU''^CiU''%U . (A.4b) 
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both S{T) and S{T') in Eq. (A. 10) are just the entropy of a rigid box of length L'i^. 
Thus the entropy-to-energy relationship is just S{E) = S'^^{L'[^; E), the entropy as a 
function of the given energy in a rigid box with length L'in. 
For such a rigid box at very low energy we find 

E^-^^u^^e-^< (A.11) 

and thus 



where oof^ — 7r/L'in = is the frequency of the lowest mode. 

Now for any given N^^ we get some tu™'^. As discussed in section 4.2, the larger 
A^in we choose, the smaller a;™'^ we will have. For fixed we can choose /? 3> l/uj^^ 
and proceed as above to get a large value of L\n- However, we are really interested in 
the continuum limit at fixed temperature. If we increase Ain with /3 fixed we will find 
that u!i (and eventually an arbitrary number of the Ua) will become smaller than 1//3. 
When this happens, the approximations of Eqs. (A.2-A.4) will no longer be good. 

However, we do not expect the entropy to decrease drastically in this limit. To 
make the entropy small would require making all the frequencies large. If the frequen- 
cies were large, the approximations we have used would again become valid. Then 
we could argue as before that the entropy should be large. It would be hard to have 
a consistent picture. 

Now consider the limit as T — > 0. For each T we start with some initial number of 
oscillators A^^^ We choose N^^^ not too large, such that uj^^ » T. With this value 
of Ain, we find L\^^ ~ n/uj^. We then let A^in — > oo and we conjecture that the 
entropy does not change much, and thus in the continuum limit S{E) ~ S^^{L'i^^; E). 
As we decrease T we can decrease the initial uj\^ and so increase L';^'^) without bound. 
Thus we make the following conjecture: 

For a given energy i?, let L'n(£') be the length of a rigid box such 
that the vacuum-bounded state with energy E and length L'^ has entropy 
S{E) = S'^{L\n{E);E). Then 

lim = oo . (A.13) 

E^O Lin 

To support this conjecture numerically we turn to direct calculation of energy 
and entropy values for vacuum-bounded states at low temperature. For various fixed 
values oi (3 = 1/T and for various numbers of oscillators we compute S and E and 
from them the equivalent length L'in. The results are plotted in Fig. A-1. While A^in 
is still small enough for the approximations Eqs. (A.2-A.4) to be valid, L'in grows 
with Ain. Once A'in has left this regime, it appears that L'in levels off. It is at least 
reasonable to believe that there is no further change in L'in as A^in — > oo. In Fig. A-2 
we plot the eventual level of L'in versus /3. It appears that the hmiting value of L'in 
grows nearly linearly with /3, and thus L'in — > oo as L? — > as conjectured. 
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Figure A-2: The length L'in of a rigid box with the same S{E) plot- 
ted against (3. Each point is the value for the largest number of 
oscillators available. 
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